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(Publication of this paper, substantially in its present form, had nearly been accomplished in July, 1941, when military con- 


siderations held it up because of beginning of actual jet propulsion in the United States. 


However, the matters presented are not 


connected with any actual development and relate only to mathematical items that might have been stated at any time since Newton.) 


INTRODUCTION 


hn OPERATION of any usual airplane requires that 
an appreciable quantity of air be taken on board 
the plane from the atmosphere and that this air later 
be discharged. The air may be used for combustion in 
the engine, for cooling of the engine or other parts, for 
a propulsive jet, or for some other purpose. Frequently, 
the intake of air is by means of a suitably located 
ram whose diffuser increases the static pressure of 
the air in a supercharger or engine inlet or provides 
pressure to force air through a cooling duct. At dis- 
charge, air or exhaust gases may be directed rearwardly 
and so exert a propelling force on the airplane. The 
limit of this is an airplane driven by propulsive jets 
rather than by a propeller. 

The intake of air exerts a retarding force on the air- 
plane, a theoretical value of which may be found by the 
application of Newton’s laws of motion. A theoretical 
value of the propelling force of a jet may be found in a 
similar fashion. These forces exist whenever air is 
taken on board the plane and air or exhaust gases are 
discharged from it, regardless of the use on the airplane. 
It is the purpose of this article to give theoretical values 
of the forces, powers, and energies for these cases and 
to discuss them briefly for an airplane flying at constant 
altitude without acceleration. 


NOMENCLATURE 


F, = 
F, = 


drag due to inlet ram, Ibs.-force 
propelling force due to jet propulsion, Ibs.-force 
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F, — F,, net thrust force 
Flow of air taken on airplane, lbs.-mass per sec. 
= Flow of exhaust gases expelled from airplane, lbs.-mass 
per sec. 
: airplane velocity, ft. per sec. 
theoretical nozzle velocity of jet with respect to the 
airplane, ft. per sec. 
dimensionless proportionality factor, equal to the stand- 
ard absolute value of the acceleration of gravity 
power required for intake, ft.lbs. per sec. 
power produced by jet propulsion, ft.lbs. per sec. 
P2 = P,. 
= energy acquired per lb. of air because of ramming in- 
take, ft.lbs. per lb.-mass 
energy to propel airplane per Ib. of exhaust gases, ft.- 
Ibs. per Ib.-mass 
E, — E, = net energy used to propel airplane, ft.lbs. 
per lb.-mass 
kinetic energy lost because of absolute velocity of ex- 
haust gases, ft.lbs. per lb.-mass 
= mechanical energy added per lb.-mass of air 
E/E», 
specific heat at constant pressure 
static temperature of the atmosphere, degrees R 
density, Ib.-mass per cu.ft. 
mechanical equivalent of heat 
flow, cu.ft. per sec. 
area, sq.ft. 
external ambient pressure 
pressure after ram diffuser 


NEWTON’S LAW RELATIONS AT INTAKE 


In order to get air on board an airplane it usually is 
necessary to accelerate it from the atmosphere up to 
the airplane velocity, a steady-state flow of air being 
produced. Consider a mass dW, of air being taken on 
board in time dt. This air must have its momentum 
changed from zero to airplane velocity. By the appli- 
cation of Newton’s second law of motion, the impulse is 
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equal to the change in momentum, or 
Fdt = (dW,/g)(V — 0) (1) 


Here F; is the force necessary to accelerate the mass of 
air. It is also the retarding force on the airplane due 
to taking the air on board. There may be substituted 
in Eq. (1) the steady-state relation w,; = dW,/dt, where- 
upon — 
F, = (w:/g)V (2) 
The power used is 
P, = KV = (wi/g)V? (3) 


The energy received by each pound-mass of air, in foot- 
pounds, is 


E, = V?/g (4) 


In the derivation of the preceding equations nothing 
whatever was said as to the location on the airplane 
where the air was to be taken on board or through what 
size or shape intake this was to be accomplished. All 
that was said was that a certain quantity of air is ac- 
celerated up to the airplane velocity. It follows that 
power of the amount given in Eq. (3) is a primary item 
regardless of how or where the air is taken on board. 
Hence, power of this amount is a reasonable datum 
against which the power used by an actual intake may 
be compared. 

These Newton law relations sidestep the details of 
the way the various pressures produce the forces, and 
this is elaborated upon in the Appendix. 


RAMMING INTAKE 


A ramming intake is an opening that faces in the 
forward direction. The purpose of this arrangement is 
to obtain a pressure beyond the intake higher than the 
static pressure of the atmosphere by making use of the 
relative motion to convert velocity into pressure. 
Since the Newton’s law relations apply regardless of 
where or how air is taken on board, it follows that they 
must apply to this case of a ramming intake. Thus, 
there is no theoretical reason why a ramming intake 
should require any more power than any other intake. 
Since an increased pressure is produced, there is every 
reason why such an intake should be used. 

Eq. (4) stated that air taken on board receives energy 
of amount V?/g ft.lbs. per lb. of air. This assumes that 
the air is brought completely to rest relative to the air- 
plane. When so brought to rest, the air has a velocity 
V with respect to the original atmosphere, so each pound 
of air has a kinetic energy V?/2g. This is just half of 
the energy, V?/g, that the air receives at intake. The 
other half, V?/2g, is potential energy that the air re- 
ceives, shown as rise in temperature and pressure. 

The mathematical fact, Eq. (4), that the air receives 
both kinetic and potential energy at intake, and re- 
ceives these in equal amounts, may be shown physi- 
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cally as well. The air just previous to intake has a 
velocity V with respect to the airplane. When the air 
is finally on the airplane, this relative velocity is prac- 
tically eliminated. Assume for the moment that this 
reduction of relative velocity takes place in a diffusing 
passage. Then, as the relative velocity is reduced, the 
pressure is increased. If the relative velocity is brought 
to a negligible value, then the air has acquired a poten- 
tial energy of amount V?/2g, and this energy will be in 
the form of increased temperature and pressure. At the 
same time that the relative velocity is being reduced, 
the absolute velocity of the air with respect to the at- 
mosphere is being increased. Thus, simultaneously 
with the acquisition of the potential energy, V?/2g, the 
air acquires a kinetic energy of the amount V?/2g. 

The form in which the potential energy of amount 

V?/2g appears needs further explanation. Assume for 
the moment that an uncompressed fluid is taken on 
board any moving vehicle. An illustration of this 
is a moving boat taking water on board. If the intake 
and diffuser has no friction or shock losses, then it is 
known from Bernoulli’s equation that the potential 
energy V?/2g would appear as a pressure rise 

pi — po = pV?/2g (5) 
when the water is brought completely to rest with re- 
spect to the boat. If there are friction and shock losses 
at intake and diffuser, all or part of the energy V?/2g is 
dissipated as heat, and the temperature of the water is 
raised. In this case the pressure rise is only a fraction 
of the value pV?/2g. In the limit there is no pressure 
rise because of diffusion, and the energy V?/2g appears 
wholly as heat. 

The intake of air onto an airplane at constant altitude 
without acceleration differs somewhat from this, since 
the air is compressible. If energy V?/2g is added toa 
perfect gas, the temperature of the gas is raised an 
amount : 


Ti oe To = V?/2gS cp (6) 


and this temperature rise takes place regardless of any 
pressure changes. If the intake and diffuser have no 
friction or shock losses, then the isentropic pressure 
rise should be produced. By equating the isentropic 
work of compression to V?/2g, the final pressure is 


pi = poll + (V2/2gIcpTo) *~* (7) 


In case there are friction and shock losses, the final 
pressure will be less than the above value. However, 
the temperature rise will be the same, since this is de- 
termined by the amount of the potential energy added 
to the gas. This is the same whether there is a perfect 
diffuser, loss of all of the velocity by shock, or anything 
between. 

For the range of airplane speeds sometimes encoun- 
tered, the pressure rise may be found approximately by 
assuming the air to be uncompressed. Thus, for an air- 
plane velocity of 400 m.p.h. at sea level, the pressure 
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rise for a frictionless and shockless intake, calculated 
for N.A.C.A. standard atmosphere for a compressed 
fluid, is 6.19 in. Hg, while the pressure rise for an un- 
compressed fluid is 5.78 in. Hg, a difference of 6.7 per 
cent. Since this difference is not great, an uncom- 
pressed fluid will be assumed throughout the rest of 
this paper. If it is desired to evaluate the effect of com- 
pressibility, Eq. (7) may be put in the form 


pi — po = cpV?/2g (8) 


where the factor c could be found by combining Eqs. 
(7) and (8). 

Since the pressure rise p V?/2g or its isentropic equiv- 
alent is produced on the average by a frictionless and 
shockless intake, it seems reasonable to use this value 
as the ideal with which the pressure rise produced by 
actual intake may be compared. This is probably the 
greatest pressure rise that can be produced by a simple 
ramming intake. It is, of course, possible to place an 
intake directly behind a propeller and in this way to 
utilize to advantage some of the irregular velocity and 
pressure disturbances created by the propeller. How- 
ever, except for an inlet influenced by the slipstream, 
from the propeller, it seems reasonable that the pressure 
rise, p V?/2g, is the greatest that can be produced with 
aram with negligible final duct velocities with respect 
to the plane. This is justified further by known pres- 
sure and velocity distributions about airfoils. It is well 
known that the total pressure about an airfoil is the 
static pressure of the atmosphere plus the velocity 
head, pV?/2g, where V is the plane velocity with respect 
to the atmosphere. Such a total pressure is often found 
everywhere around an airfoil except on the top side near 
the rear. Here a lower pressure may be found, since 
the flow is undoubtedly influenced by the boundary 
layer. Thus, no case of a static pressure rise greater 
than pV?/2g is known to be produced by simple ram. 

It has sometimes been stated that the forward stag- 
nation point at the leading edge of the wing section is 
the optimum location for a ramming intake, since at 
this point the velocity head has been converted into 
pressure. This assumption has gone so far as to sup- 
pose that there could be flow with velocity V into an in- 
take, starting with the static pressure that would have 
existed when there was no flowinto the intake. This is 
greater than the ambient pressure by pV?/2g. Then, 
after a perfect ram, there would be a static pressure of 
pV?/g. However, the stagnation point is the place on 
the forward edge of the wing where the flow divides to 
pass on either side. There is no flow at this point, and 
the static pressure is increased by the amount pV?/2g. 
This is identical with the action of an impact tube. At 
an intake there is flow of the air, so that, if an intake is 
located where the stagnation point would otherwise be, 
it follows that the stagnation point must change its 
location. Therefore, it is a question as to whether or 
not the stagnation point as a location for an intake has 


any special advantage over any other point on a wing. 
However, there is need for further ram research. 

It has been shown that a pressure rise of pV*/2g is 
theoretically possible, but that this may be reduced by 
friction or shock. These factors depend upon the flow 
area and shape of the intake, as well as upon the loca- 
tion. If the intake is at the nose of a wing section or 
fuselage or at.a forward scoop projecting some distance 
from the airplane, the area necessary to get a flow Q on 
board is 


A=Q/V (9) 


where V is the airplane velocity. If the air is to be used 
for combustion, the flow Q is determined principally 
by the motor characteristics, so that in this case the 
flow Q is fixed chiefly by conditions beyond the intake. 
By means of Eq. (9) an area may be found so that the 
air has a relative velocity V with respect to the airplane 
at intake. Since this relative velocity is to be reduced 
and a pressure rise produced, in accordance with Ber- 
noulli’s equation, the passageway should be, theoreti- 
cally, a diffusing one. However, it is a matter for ex- 
periment to determine whether the area given by Eq. 
(9) is best and whether a diffusing passageway does 
produce the maximum pressure rise. Flow conditions 
at an intake are such that some other area or form may 
be found to be preferable, although in no case would a 
pressure rise greater than pV?/2g be expected. Further- 
more, with the actual flow Q fixed by an engine displace- 
ment or by nozzle area of a jet propulsion gas turbine, 
the area required for a given V and Q as given by Eq. 
(9) will not, in general, match the actual area of the 
ramming intake. If the actual area is too great, some 
of the flow heading for the ramming intake will not 
enter it but will branch off at the edges. If the ram- 
ming intake area is too small for a given V and Q, 
other losses will occur which are not here investigated, 
these probably needing research. 

Consideration must be given not only to pressure rise 
but also to retarding forces on the plane. It has been 
shown that the Newton’s law force at intake is w,V/g. 
It seems reasonable to expect that any intake will 
affect the flow of air around a plane and change the 
drag. The amount of this change depends on the 
location and form of the intake. It is not apparent 
whether the drag is increased or reduced, and it is im- 
practical to try to evaluate the change theoretically. 
The Newton’s law retarding force on an airplane be- 
cause of an intake is w,V/g, but the total effect of this 
retarding force, including the change in the airplane 
drag, remains undetermined. It is a matter for experi- 
ment to determine the optimum location and form for 
a ramming intake so that the total force holding the air- 
plane back is a minimum and so that as large a pressure 
rise as possible is ohtained. A similar situation exists at 
discharge of cooling air, exhaust, or a propulsive jet. 
The Newton’s law forces for intake and discharge are 
imposed on an airplane, which, without the openings 
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for the intake and discharge, has a certain drag. The 
changes in the airplane drag due to the presence of the 
openings may be called “intake interference drag’’ and 
“discharge interference drag.” 

Finally, any intake that supplies air for combustion 
should be in a region of steady flow. If the intake 
should be in a region of pulsating flow and there is a 
compressor beyond, the pulsations may, be magnified 
and produce an undesirable effect. 


NEWTON’sS LAW RELATIONS AT DISCHARGE 


The theoretical forces produced by a rearward dis- 
charge of air or exhaust gases from an airplane may be 
found in a manner similar to that with which the forces 
at intake were found. Suppose that a mass of gases w» 
is discharged rearwardly with a velocity relative to the 
airplane of V2. These gases can have no knowledge 
as to whether the airplane is at rest or in a condition of 
steady motion. They merely have their velocity rela- 
tive to the airplane changed from a negligible value up 
to the velocity V2. Thus, as in Eq. (2), the force acting 
on the gases to accelerate them from an original neg- 
ligible velocity with respect to the plane up to the veloc- 
ity Vo is 

F, = We V2/g (10) 


If the gases are discharged rearwardly this force acts 
on the airplane in a forward direction and, hence, helps 
to propel it. The propulsive power used in this way is 


P, = F,V = w.V2V/g (11) 


The propulsive energy from each pound-mass of gas 
discharged is 


PROPULSIVE JETS 


Thus a rearward discharge exerts a propulsive force 
on an airplane. Such a discharge may be the engine 
exhaust gases, the engine cooling air, a propulsion jet, 
or a discharge for any other reason. The theoretical 
forces, powers, and energies arising from consideration 
of only the discharge have been given. It should be 
noted that the propulsive force involves the jet speed 
relative to the airplane but does not involve the air- 
plane speed. For this reason a propulsion jet is some- 
times considered as a constant thrust device. 

A propulsive jet undoubtedly affects the parasitic 
drag on the airplane, as is the case of ramming intake, 
but the amount of this change cannot be predicted by 
simple theory. It has been suggested that the jet 
could be placed so that it would help clear away the 
boundary layer on the top rear part of a wing, and in 
this way improve the performance of the wing, or that 
a low pressure region could be found so that the jet 
velocity would be greater than that obtained by an ex- 
pansion to the static pressure of the atmosphere. It is 
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a matter for experiment to determine whether there is 
any way to discharge a jet with a velocity V2, deter- 
mined by expansion to atmospheric pressure, without 
extra drag due to the jet or with no extra drag and a 
greater velocity. Until such knowledge is available it 
seems best to assume the theoretical values given 
above. 

The fact that the propulsive energy per pound-mass 
of gas discharged is given by the expression V2.V/g may 
be given a physical interpretation by means of the 
mathematical identity 

2 7,2 = 2 
Ee eee eee 
g +g g “8 

The quantity V?/2g is the kinetic energy of the gases 
while they are yet on the airplane in a duct with negli- 
gible velocity. In order to discharge them from the 
airplane it is necessary to give them a velocity relative 
to it, which is obtained by a conversion of a pressure 
produced on the airplane, into relative velocity. That 
is, previous to discharge, the gases must have a pressure 
greater than the static pressure of the atmosphere. 
This pressure represents potential energy V2?/2g, con- 
verted into kinetic energy by the jet nozzle, and the 
gases have this in addition to the kinetic energy, V?/2g. 
These two items are directly comparable to conditions 
at an intake where the air acquires both kinetic energy 
and pressure. However, at discharge a third item is 
involved. The relative velocity of the gases at dis- 
charge generally differs from the airplane velocity, so 
that*the jet has an absolute velocity in space, (V2 — 
V), with respect to the atmosphere. Since this velocity 
is eventually dissipated by friction, it represents kinetic 
energy (V2 — V)*/2g, which is entirely wasted. To ob- 
tain the energy used to propel the airplane it is neces- 
sary to subtract this loss from the kinetic and potential 
energies the gases have just previous to discharge, as is 
done in Eq. (13). 


PROPULSIVE JET WITH INTAKE FORCE DEDUCTION 


Consideration has been given to the forces, powers, 
and energies of discharge. If the gases were originally 
on board the airplane, as in the case of a rocket, this 
would be the complete procedure. However, often, 
air has to be taken on board the airplane from the 
atmosphere. To find a true measure of the effectiveness 
of a propulsive jet, it is necessary to subtract the retard- 
ing force and the power at intake from the jet propulsive 
force and power at discharge. Thus the net propulsive 
force is 


F = F, ial F, = (w/z) [ (we Vo/w1) = V) (14) 
This equation completely sidesteps the pressures on the 
various walls that serve to produce F, and these are de- 
duced in detail in the appendix. The net power used to 
drive the airplane is 


P = P, — Py = (wiV/g)[(w2V2/w:) — V) (15) 
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The net energy used by the airplane per pound-mass of 
air taken on board is 


E = E, — E, = (V/g)[(weV2/uw1) — V] (16) 


These equations apply to any backward discharge from 
an airplane. In the case of a radiator or cooling duct, 
the discharge flow would be equal to the intake flow, 
W2 = w,, and the above equations are somewhat sim- 
plified. In the case of air taken on board to be used for 
combustion in a power plant, the discharge flow would 
be greater than the intake by the amount of fuel added. 
For the case of a nozzle discharging a jet of exhaust 
gases, it may be assumed that V2 is an ideal velocity 
and that the actual relative velocity is less than this by 
a nozzle velocity coefficient or other loss factor w/w», 
which is a reasonable value. Then, for discharge of 
cooling air, exhaust, or a propulsive jet, Eqs. (14) and 
(16) will be written 


F= F, — Ff, 


(w/g)(V2 — V) (17) 
and 
E=f-£ (V/g)(V2 — V) (18) 


It is interesting to note that the net propulsive force 
and energy given by Eqs. (17) and (18) depend on the 
velocity difference (V2 — V). This is the final absolute 
velocity of the gases. If V2 = V, the jet is stationary 
with respect to the atmosphere so that the propulsive 
force of the jet is equal to the retarding force due to 
intake. If the propulsive force is to be greater than the 
force at intake, it follows that the jet velocity relative 
to the airplane must be greater than the airplane veloc- 
ity. This produces an absolute velocity giving the kine- 
tic energy loss 


E; = (V2 — V)?/2g (19) 


Thus, for a jet to produce a propulsive force at discharge 
greater than the retarding force at intake, there is an 
inevitable loss in the kinetic energy of the jet. It 
should be noted that the relative velocity at discharge 
may be less than the airplane velocity. In this case 
there is a propulsive force at the discharge, given by 
Eq. (10), but the retarding force at intake is greater 
than the jet force at discharge so that the net effect on 
the airplane is a retarding force as given in Eq. (17). 

To produce a jet an energy input is necessary, which 
is given by 


E, = (1/2g)(V2? — V*) (20) 


This quantity E, is the potential energy added to the 
gases while on the airplane. It has already been shown 
that previous to jet discharge the gases on the airplane 
are at a pressure above that of the atmosphere and 
that this pressure energy is given by V2?/2g. Likewise, 
in the discussion of ramming intake it was shown that 
the air acquired pressure energy, V?/2g, by means of 
the ramming intake. Hence, the difference (V2? — 
V*)/2g is the potential energy added to the air while it 


is on the airplane. Nothing is said here as to how this 
energy is added. 

There has now been discussed the propulsive energy 
E, the kinetic energy loss E3, and the net potential 
energy added to the gases on the airplane, E,. The 
kinetic energy loss is the only loss. Hence 


E+E;=E, (21) 


This relation may be verified by the substitution of 
Eqs. (18), (19), and (20). Since the input is Z, and 
the output is the j. opulsive energy £, the ideal effi- 
ciency is given as 


n = E/E, = (2V/V2)/[1 + (V/V2)] 22) 


It can be seen that » is a function of the ratio V/ V2 and 
is so plotted in Fig. 1. This shows the reduction in effi- 
ciency when the airplane speed is less than the jet 
speed, although this condition is necessary if the pro- 
pulsive force at discharge is to be greater than the re- 
tarding force at intake. 

So far, the jet speed has been related to the airplane 
speed, but it remains to find the best value of the air- 
plane speed. To do this it will be assumed that a con- 
stant amount of potential energy is added to each pound 
of air as it goes through the airplane. Thus, in Eq. 
(20) it is assumed that E, is held constant. Then 


V. = ¥2eF, + V? (23) 


and the efficiency 7 can be found in terms of the air- 
plane speed. 


n = 2V[¥2eE, + V? — V]/2gE, (24) - 


This equation could be plotted, but this is not here in- 
cluded. The best airplane speed can now be found by 
setting the derivative dn/dV to zero. Thus 


dy _ [Vez + - vi) 
dV gEy\2gE, + i 








0 (25) 


or 
29, + V2 = V2 (26) 


This equation can be satisfied only if E, is zero or V is 
infinite. If E, is zero, there is no net propulsive effect ’ 
on the airplane. It follows that, theoretically, the best 
speed for an airplane with jet propulsion is infinite. 
Such a speed is, of course, impossible, but the analysis 
shows that for good propulsive efficiency the airplane 
speed must be high. 
The thrust that is obtained by jet propulsion is given 

in Eq. (17), 


F = (w/g)(V2 — V) (17) 
If Eq. (20) is substituted in this, there results 
F = (w/g) [ y2e2> +Ve— v| (27) 


The thrust when the airplane is at rest is found to be 
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Propulsion efficiency and thrust of a jet. 
= 2(V/V;) 
1 + (V/Vs2) 
gh V 1 — (V/V2) 
Fo 1 + (V/V2) 


plane velocity 


Fic. 1. 


V 


V2 = jet velocity relative to plane 
Fy = jet thrust at take-off 
F = jet thrust at velocity, V 


In deriving the expression for F/ Fy it is assumed that a constant 
amount of energy Ep = (1/2g) (V2? — V?) is added to each pound 
of air passing through. 


Fy = (w/g) ¥2E, (28) 


Note that the quantity ~/2gE, is the jet velocity when 
the airplane is at rest. It is convenient to use this 
quantity when considering the thrust produced when 
the airplane is moving. Thus 


F ; v2 V 1 — (V/V») 
cae eA tb Cian i ee Ete OBE 
Fo \ * 2gE, V2cE, ' + (V/Vs) ia 
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This equation is plotted in Fig. 1. It can be seen that 
the thrust is a maximum when the airplane is at rest— 
in other words, at take-off. This is a fortunate cir- 
cumstance, since it is at take-off that the maximum 
thrust is needed. This occurs if the airplane takes on 
air as rapidly at take-off as it does later. This is prob- 
ably not the case, but the thrust produced for each 
pound of air going through the airplane is a maximum 
at take-off. 

Before concluding the discussion of jet propulsion, it 
is in order to compare the performance of jet propul- 
sion with that of propeller propulsion. It is well known 
that to produce a thrust a propeller must create a stream 
of air going backward with a velocity greater than the 
forward speed of the airplane, thus giving rise to a ki- 
netic energy loss. This loss is similar to the kinetic 
energy loss with jet propulsion. There is, however, 
one difference between the two. A propeller generally 
stirs up a large amount of air and adds but little energy 
toeach pound of it. Hence it is suitable for low velocity 
travel. A propulsive jet of high velocity can be pro- 
duced, and such a jet is suitable for high velocity travel. 

Finally, it should be noted that no consideration has 
been given to the effect-of a propeller slipstream on an 
intake and discharge, such as occurs in the N.A.C.A. 
cowling or other openings. 

The authors gratefully acknowledge the assistance 
given in writing this paper by Dale D. Streid, Charles T. 
Oergel, and Dr. Chester W. Smith. 
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Appenpix— Pressures ON W atts OF A Vesset Whurcr Provive THE Jer 
Momentum Forces 


(Added in 1945 with the collaboration of C. W. Smith) 


Consider a vessel moving without acceleration, or 
stationary. This may be a jet-propelled airplane at 
constant altitude, a chamber for measuring nozzle re- 
action, a bend in a pipe carrying flowing fluid, or any 
similar vessel. An axial direction is defined, for con- 
venience, parallel and opposite in direction to a principal 
jet leaving the vessel through an opening. There is 
supposed a variable internal pressure p, and an external 
ambient pressure fp on all external walls, and the dis- 
charge of a fluid jet with the same final pressure pp, 
through an area A» with a velocity V2. There may be 
any number of other inlet and outlet openings, and the 
area A to be considered is the total projected area of all 


fluid openings in planes perpendicular to the axial di- 
rection. 

A major use of the equations herein is the measure- 
ment of nozzle velocity by weighing nozzle reaction. 
In this case the fluid inlet area is in a plane parallel to 
the axial direction and so is taken as zero. Another use 
is computation of force on a pipe elbow; another is ef- 
fect of air taken on and discharged from an airplane 
flying without acceleration at constant altitude (per- 
haps with propulsion jet and an augmenter) and with 
V2 less than sonic velocity or with a nozzle expanding 
portion suited to supersonic V2. The general purpose 
of this appendix is to show the parts played by the fluid 





pre 
th 
ba 


40 

not 
Fig 
anc¢ 


RELATI 
INLET 


Fm 


F; 


3 


el ie 


f(V) 


The 
numbe 
po, So t 
be res¢ 
at any 
becaus: 
an axia 
gives a 
force is 

Insta 
indeper 
the bod 
interna 
and box 
b, whic! 


hat 


cir- 
um 

on 
ob- 
ach 
um 


1, it 
yul- 


am 
the 


etic 
yer, 
ally 
Tgy 
ity 
yro- 
vel. 
has 
| an 
ee 


nce 
pe 


; On 
431, 


take 








RAMMING INTAKE AND PROPULSIVE JETS 117 


pressures on the vessel walls. These are sidestepped in 
the previously given and usual treatment of the problem 
based on theory of conservation of momentum of the 
flowing fluid. (The senior author remarks that during 
40 years of work with nozzle reaction, etc., there has 
not been the clear picture shown by the analysis herein.) 
Fig. 2 gives a diagram for the particular case of inlet 
and outlet flow parallel to the axis. 





EXTERNAL PRESSURE Po ON ALL WALLS, GIVING 
A NET EXTERNAL AXIAL FORCE, Po(A)~A>) 
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es | Az —» V2 
re: a a yates =I RELATIVE 
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—- dA PROJECTED AREA 





——F NET AXIAL THRUST 
ON ALL VESSEL 
WALLS DUE TO 
INSIDE AND OUTSIDE 
PRESSURES 


Fm = FORCE ACTING TO CHANGE THE MOMENTUM 
OF THE FLUID WITHIN THE. VESSEL AT 
ANY INSTANT 

—— = POSITIVE AXIAL 
DIRECTION 


Fic. 2. Diagram in a plane parallel to the axial direction. 


NOMENCLATURE FOR APPENDIX 


JS pdA = internal pressure axial force, equal and opposite 
for vessel and fluid 

w = flow (lb.-mass per sec.) 

Po = external ambient pressure (constant) 

p = internal pressure (variable) 

A, = sum of all inlet areas projected on planes perpen- 
dicular to the axis 

Az = exit areas projected on planes perpendicular to the 
axis 

A = projected area in general, in a plane perpendicular 


to axial direction 

= inlet velocity 

= exit velocity 

= net force acting on fluid 

= net force acting on walls 

= dimensionless factor equal to acceleration of 
gravity 

(VY) = drag on outer walls (which do not include A; and 

A») due to vessel motion 


sd a ta 


The vessel walls inside and out are equivalent to a 
number of steps on each of which is the pressure p or 
po, so that the axial force is p dA or podA. There could 
be resolved, in the axial direction, pressures on a wall 
at any angle, which would give the same result. This is 
because a pressure fp on an inclined wall area ds gives 
an axial force po ds cos 6. But an inclined wall area ds 
gives a projected area dA = ds cos 6. Hence the axial 
force is po dA or pdA. 

Instantaneous force equations will be written for two 
independent bodies: first, the vessel itself; second, 
the body of fluid within the vessel. There is a common 
internal force acting equally and oppositely on vessel 
and body of fluid due to the variable internal pressure 
p, which is given by fp dA. 


In a jet-propelled plane the vessel that is considered 
is the passage system through which the fluid passes 
from inlet to outlet. It has various wall areas on which 
this pressure p is exerted. In the case of a plane with 
ramming intake and diffuser, there is a positive force 
on the diffuser walls. In a jet-propelled plane the pas- 
sage walls include the air compressor chamber and the 
combustion chamber, where there may be compara- 
tively high pressures. Finally, there is the jet discharge 
nozzle within which the pressure gradually decreases to 
the ambient value ~p. The net difference between the 
high pressure on the inner wall area and the gradually 
decreasing pressure on the nozzle walls is the major 
force propelling a jet plane. This is visualized most 
easily by thinking of a vessel used for nozzle reaction. 
Here the nozzle hole is a gap in the wall, so there is a 
thrust force due to the corresponding area at the oppo- 
site side of the vessel, as well as decreased pressure near 
the discharge nozzle. (All of the internal space in the 
plane, other than the complete passage between inlet 
and discharge jets, constitutes a completely closed ves- 
sel whose interior is filled with air at the ambient pres- 
sure fp. There is no net axial force from the pressure 
within such a closed vessel.) 

In addition to the equal and opposite interior pressure 
on vessel and body of fluid, there are exterior forces due 
to the openings in the vessel walls. For a vessel for 
measurement of nozzle reaction, or a 90° elbow, there is 
one opening in a plane perpendicular to the axial direc- 
tion, where the jet is discharged. The inlet opening 
is zero in a plane perpendicular to the axial direction. 
For an airplane with ram inlet, or a jet-propelled plane, 
there are inlet and outlet openings in planes perpen- 
dicular to the axial direction. If there were no such 
openings in the vessel, the net force on the outer walls 
due to the ambient pressure fp would be zero. For a jet 
propulsive apparatus on block test, there is a large inlet 
opening and a smaller opening for the jet discharge, as 
shown diagrammatically in Fig. 2, so that there is a net 
external force pushing the vessel in the axial direction, 


po(A1 — Az) 


This is a real part of the force on the vessel, whose ex- 
istence often is forgotten. For the case of a vessel for 
measuring nozzle reaction, A; is zero, and the force is 
—poAz. This is the force on the vessel due to the pres- 
sure on_ all of the external walls but with no wall at the 
nozzle opening. For the case of an airplane with cooling 
air entering and leaving, A, may be larger than A; be- 
cause of the lower velocity, and so f(A: — Az) may bea 
retarding force. In any case the effective vessel surface 
on which the external pressure acts is the difference 
between the entrance and exit areas A, and A. Itis to 
be noted that we assume the areas A; and A; to be com- 
pletely filled by continuous flow of fluid, with the veloci- 
ties V and V2. Otherwise complications arise not here 
investigated. No method is here proposed for accom- 
plishing this exact filling at all speeds and all altitudes. 
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For a jet propulsive apparatus on block test the am- 
bient pressure gives the above axial force on the differ- 
ence of inlet and jet openings. When the jet-propelled 
plane is in motion, there is a retarding drag force f(V). 
This is defined as being the resistance introduced by 


motion. 
So much for the pressure forces on the external vessel 


walls due to the inlet and outlet openings. There also 
are forces on the fluid passing through the vessel due to 
the ambient pressure of the inlet and outlet fluid 


streams 
—poAi + pode 


or 
—po(Ai — Az) 

This is exactly equal and opposite to the force on the 
external vessel walls due to the inlet and outlet open- 
ings. 

So the gross forces F and F,, on the vessel walls and 
on the fluid within the vessel have the numerically 
equal but opposite terms, with absolute magnitudes of 


JS p dA and p(A; — A2) 
that is, 
F= fpdA + po(Ai — A:) 
and 
F, = —F 


The force F on the vessel walls is the real force that 
pushes it and should have further analysis. One part, 
po(A1 — Ag), is the ambient pressure, po, acting on that 
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part of the external walls not neutralized by an equal 
wall area on the opposite side because of an opening. 

Inside of most of the vessels being considered, there 
is for the most part a high pressure p, except near the 
nozzle opening. Here the pressure on the nozzle walls 
is gradually reduced to the ambient value while fluid is 
being accelerated. The pressure p on the corresponding 
opposite wall gives an internal pushing force. The total 
force is the sum expressed mathematically by the value 
of F. This is written most easily by computing the 
numerically equal value of F,, from the Newton’s law 
forces as already done, giving 


F = (w/g)(V2 — V) 


This sidesteps the details in this appendix but misses 
the analysis given. In the case of a jet-propelled plane 
without acceleration, at constant altitude, F = f(V), 
where f(V) is the drag as defined herein. 

An interesting detail is the force of nozzle reaction 
with an uncompressed fluid, which is computed from 
the jet momentum and is given in many textbooks as 


F = 2A(Ap) 


That is, the force is the pressure difference on twice the 
unbalanced area, because of the decreased internal pres- 
sure on the internal walls near the nozzle opening. It 
is to be repeated that in all of the matter in this paper 
relating to a jet-propelled plane, we assume constant 
altitude and no acceleration and inlet and outlet open- 
ings are assumed to be filled with fluid at the ambient 
pressures pp. This postulates exactly the proper areas 
and a proper nozzle, with expanding portion if there is 
supersonic velocity. 
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Stability of Semimonocoque Wing Structures’ 


HENRY L. LANGHAARt 
Consolidated Vultee Aircraft Corporation 


ABSTRACT 


The Consolidated Vultee Aircraft Corporation recently built 
and tested several cambered box beams that were structural mod- 
els of stressed-skin airplane wings. It was found, as a result of 
these tests, that plate-stringer column-fixity factors of the speci- 
mens ranged from 1.00 to 3.00. These observations are rather 
revealing in view of the frequent arguments among stress analysts 
concerning the correct fixity factor to be used for calculating 
allowable plate-stringer loads. It is therefore especially satisfy- 
ing to find that the test data are in harmony with a simple ap- 
proximate theory of elastic stability of wing structures. This 
theory exposes the possibility of occasionally exerting control over 
fixity factors for purposes of weight economy or bulkhead re- 
arrangement. It further suggests modifications of wing con- 
struction, such as the substitution of shallow ribs for bulkheads 
in highly cambered wings. 


SYMBOLS 


= bulkhead spacing 

peripheral interspar distance 

= skin thickness 

te = area of a representative unit width of the plate-stringer 
cross section. (The area ¢, may be called ‘‘equivalent 
thickness’”’ because it is the thickness of a homogene- 
ous shell with the same weight and area as the given 
plate-stringer panel) 

p = radius of gyration of an element of the plate-stringer 
cross section with reference to its centroidal axis 

Z = ordinate of the centerline of the plate-stringer cross 

section measured from the neutral plane of the wing 

mez, = the maximum value of z 

= a mean value of z, defined by Eq. (14) 

peripheral and spanwise coordinates (see Fig. 1) 

ax/B,n = wy/L, ¢ = 2/B (dimensionless coordinates) 

= radius of curvature of a cross section of the compression 

panel at the point x 

= B/xR, a dimensionless measure of the curvature 

the panel fixity constant 

m =the number of half waves in a cross section of the 
buckled panel 

r = the ratio of the elastic modulus of a bulkhead chord 
member to the effective modulus of the plate-stringer 
combination 

y = Poisson’s ratio 

J = torsional rigidity of a bulkhead chord member 

a,‘™ = coefficients in the Fourier analysis, Eq. (21) 

A, A;, Ao, K = dimensionless constants defined by Eq. (10) 
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INTRODUCTION 


pon COMPRESSION SURFACE of a semimonocoque 
wing is structurally analogous to a set of parallel 
columns that acquire compression loads by the bending 
of the wing. Accordingly, a column fixity factor C is 
commonly employed for determining allowable stringer 
stresses. There is, however, a tendency by some engi- 
neers to discredit the use of a fixity coefficient C for 
short-column calculations—a bias that seems to stem 
from the idea that the reduced column length L’ has a 
more fundamental significance than the constant C. 
This, of course, is not the case if one adheres to the 
usual definition L’ = L/+/C or to the general Euler 
column formula with a reduced modulus. The main 
reason for preferring C rather than L’ in the present 
study is that the former constant is dimensionless. 

There are several effects that cause the constant C to 
differ from unity. First, there is the possibility that the 
bulkheads are so light that they participate in the buck- 
ling. This condition inherently contributes a negative 
increment to C. However, for the present study, the 
phenomenon of linear bulkhead deflection is excluded. 
There then appear to be only two effects that can in- 
fluence stability. These are the restraints furnished by 
the skin and by the bulkhead connections. 

The restraining effect of the skin consists essentially 
of two parts, membrane action and flexural action. The 
former action is known to be peculiar to curved shells; 
it consequently has significance only for cambered 
wings. Likewise, the flexural stiffness of the skin and 
the torsional stiffness of the bulkhead chord members 
are ordinarily negligible for flat compression panels. 
However, these factors should not be disregarded in the 
theory of stability of a cambered wing, for without 
them the mode of buckling becomes indeterminate. 
This circumstance arises because the flexural stiffness 
of the skin and the torsional stiffness of the bulkhead 
chord members do have significance for modes of buck- 
ling which present several waves in a cross section of a 
panel. 

In order to formulate the stability problem that is 
suggested by the above remarks, the view is adopted 
that the compression surface of a wing is an orthotropic 
shell. The significant plate-stringer dimensions are 
then the skin thickness /, the radius of gyration p, and 
the thickness ¢, of an equivalent homogeneous shell. 
Attention is confined to a single, cylindrical, compres- 
sion panel that is constrained so that the normal and 
transverse components of the elastic displacement vec- 
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tor vanish at the ends (i.e., at the bulkheads), while only 
the normal component vanishes at the longitudinal 
edges (i.e., at the wing spars). Rotational restraints 
of uniform bulkhead chord members are considered to 
act at the ends of the panel. 


POTENTIAL ENERGY OF A BUCKLED PANEL 


The cylindrical compression panel under considera- 
tion is schematically shown in Fig. 1. The type of edge 
support that is considered has been discussed in the pre- 
ceding paragraph. Peripheral and axial coordinates 
(x, y) are set up as shown in the figure. The bulkhead 
spacing L and the peripheral interspar distance B are 
indicated by the bulkhead locations y = 0, y = L and 
the spar locations x = 0, x = B. 








\, 


Fic. 1. Diagram of compression panel. 


The axial load P is considered to be proportional to 
the ordinate z, measured from a fixed plane. The maxi- 
mum value of P may be represented in the form 


Psa. = Cr°Et.p?/L* (1) 


in which /,p? is the moment of inertia of a unit width of 
the plate-stringer cross section and £ is Young’s modu- 
lus. Since this relation formally agrees with the 
Euler column formula, the constant C denotes the 
fixity factor that is conventionally used when margins 
of safety are based upon the most highly stressed 
stringer. 

The physical idea behind the following mathematical 
investigation is the energy theory of elastic stability, 
which characterizes the buckling point of an elastic 
system by an equality between the work of the external 
forces and the increment of strain energy in a small 
buckling deformation. A virtual deformation is de- 
fined by the transverse, axial, and normal components 
(u, v, w) of a vector displacement function (Fig. 1). It 
is assumed in the sequel that the elastic restoring forces 
exerted by the stringers preclude axial displacements v 
other than the second-order displacements attending 
bending of the stringers. With this assumption, the 


1946 


various components of potential energy are evalu- 
ated. 
(A) Virtual Work of the External Forces 
By geometry, the axial contraction of a stringer 
is 
= /2 fo" (Ow/dy)* dy 
and, consequently, the virtual work of the compression 
load is 
W = 1/2 fo? Pdx fol (Qw/dy)? dy (2) 
It is convenient to employ, instead of the coordinates 
(x, y), the dimensionless coordinates 
& = rx/B, n = ry/L (3) 
With these variables, Eqs. (1) and (2) yield 


2 2 . wv ° 
W = Cr BEt,p f F z 
2L? YS ee 


in which the subscript 7 denotes a partial derivative. 


w,” dé dn (4) 





(B) Energy of the Stringers 


The strain energy of bending of the stringers, as de- 
termined by elementary beam theory, is 


U; = /2Et.p? fo? fo'w,,? dx dy 


With the introduction of the variables (£, 7), this ex- 
pression becomes 


“. ° 2. eee fe 
U, = nc Ot RON dé dn (o) 


(C) Energy of Flexure of the Skin 


The strain energy of bending of the skin is developed 
in the general theory of elastic plates.' Specifically, 


TSE aS af a 
" 24(1 — vy?) 
[Wrr? + 2vwWp,Wyy + 2(1 — v) wry?) dx dy 


The term w,,’ is here discarded because it is contained in 
the expression for the energy of the stringers (Eq. 5). 
In terms of the variables (£, 7), the above equation be- 


comes 
B 
WerWay + 
20 Fe 


to gl Solow 
' aa — »?)B3 
2(1 — ») (7) Ws | dé dn (6) 


(D) Membrane Energy of the Skin 


Since v = 0, the membrane strain components are 
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section at the point & The fundamental membrane- 
energy expression, 


[en? + €,? + Qverey + '/2(1 — v)y?] dx dy 


then reduces to 


ss EtL 
a 2(1 arth fs  (« i = = 
1/ — yp B* 2 
Al — ») (2) wu: | ae an (7) 


(E) Energy of the Bulkhead Chord Members 





In general, the restraining effect of a bulkhead is di- 
vided between two consecutive compression panels, and 
therefore the strain energy of only one bulkhead chord 
should be assigned to a panel. The twist per 
unit length of the bulkhead chord at y = 0 is 


Wry (x, 0). 
Hence, if J denotes the torsional stiffness of the 
chord member, the strain energy is 


Us = VeGI So? wey%(x, 0) dx 
or 


Us = (x°GJ/2BL*) fo"w;,7(é, 0) dé (8) 


(F) General Equation of Buckling 


With Eqs. (4) to (8), the buckling criterion (W = 
Ui > U2 + U3 + U4) yields 


cf rh * w,2dtdn = Lf Wy? dé dn + 
a0, Suen o Jo 
af Fé [(ug + xw)? + Ku,?] dé dy + 
o J 
Af 4 (we? + WW, + 4Ku,, :) di dn + 
0 Jo -y 








~ Ag f ¢y°(E, 0) dé (9) 


where, for conciseness, the following dimensionless 


nomenclature is employed 





A = L4t/x*B’p*t,(1 — v?) 

A, = 7°A??/12B? 

A, = LJ/B*p*t,(1 + v) 

\ = ratio of Young’s modulus for a bulkhead - (10) 
chord to the effective modulus of the 
plate-stringer combination. 

« = B/(rR) 

K = 0.5 (1 — v)B?/L? J 


It may be noted that Young’s modulus cancels out of 
Eq. (9), except for the factor A, and it is therefore to be 
expected that this equation remains valid for the short- 
column range, provided that \ is calculated with a re- 
duced modulus. 


The forced boundary conditions that attend Eq. (9) 
are 


w(t, 0) = w(é, r) = w(0, n) = w(x, n) = 0 ary 
u(é, 0) = u(é, r) = 0 
In addition, there are the natural boundary conditions 


u,(0, n) = u; (7, n) = 0 (12) 


which arise from the assumption that the spar flanges 
have negligible lateral stiffness. These relations are 
simply derived from the stress-strain relation, 


= [E/(1 — v)](e + ve.) = 0 


With the initial assumption, v = 0, this reduces to «, = 
0. In view of Eq. (11), the vanishing of €, is equivalent 
to Eq. (12). 

The effect of lateral stiffness of the spar flanges may 
be incorporated in the theory as a perturbation by add- 
ing the strain energy of spar bending to the other com- 
ponents of potential energy while retaining Eq. (12). 
However, this refinement is usually insignificant. 


APPROXIMATE MATHEMATICAL THEORY 


The minimum value of the constant C that is consis- 
tent with Eqs. (9), (11), and (12) will be sought with 
the aid of the assumption that, in the process of buck- 
ling, the panel acquires a sinusoidal wave form with a 
single longitudinal loop but several chordwise loops. 
Mathematically, this assumption is expressed 


W = Wo sin mé sin 7 (13) 


A thorough discussion of the validity of this assumption 
cannot be given. However, it is pertinent to remark 
that, if the skin stiffness and the rotational end re- 
straints are disregarded and if the load distribution P 
is uniform, then the problem admits a rigorous treat- 
ment without the use of assumption (13). This leads 
to the conclusion that the fixity constant C is an eigen- 
value of the Sturm-Liouville differential equation 


¢” — Ko + [AKx*/(C — 1)]¢ = 


with the boundary conditions ¢(0) = ¢(7) = 0. The 
essential feature of the secondary restraints is that they 
furnish criteria for selecting a particular one of these 
eigenvalues. The corresponding eigenfunction ¢ is pro- 
portional to the normal deflection function w with a fac- 
tor that depends upon x. In the case of circular camber 
(x = const.) the ¢ functions are seen to be circular 
functions, and it follows that Eq. (13) is theoretically 
correct for this case. 

With the introduction of Eq. (13), the left side of Eq. 


(9) may be evaluated as follows: 
0” ", 
f z sin? mé dé 
0 


LF 
0 0 <2mar, 


Now, by the mean-value theorem for integrals, there 








w,? dé dyn = 


Zmar. 
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exists a constant 2 in the range 2min, < 2 < 2maz,, Such 
that 


= (2/r) fo"z sin? mé dé (14) 


hewn © 


The other winagaie's in Eq. (9) which depend only 
upon w may be immediately evaluated. Thus, Eq. (9) 
is reduced to 


dt)-14 
Zmaz. 

4” [Wea + sa)? + Ku!) de dn + 
T“Wo 0 0 


Am‘ + (* ) me + \Aym*? (16) 
Vv 


Hence, 

















It is now not difficult to choose the function u(é, 7) 
to minimize the integral 


Q= fT Solu; + ew)? + Ku,?] dé'dn 


For this integral, Euler’s equation of the calculus of 
variations is 


(17) 


Use + Ku,, = —(0/0€)(xw) (18) 


A solution of Eq. (18) which satisfies the boundary con- 
ditions (11) and (12) is representable in the form 


. w= f—) sing 


This form automatically satisfies the boundary condi- 
tions (11), and it reduces Eq. (18) to the ordinary dif- 
ferential equation 


f" — Kf = —wo(d/dé) (x sin mé) 
Since the function «x sin mé is determined by the 


geometry of the airfoil, it may be represented by a 
Fourier analysis of the type 


(19) 


(20) 


«sin m~ = >> a,™ sin né (21) 
n=1 
Eq. (20) then becomes 
f" — Kf = —w > na,™ cos né (22) 
n=1 


By Eq. (12), the boundary conditions associated with 
this equation are 


f'0) = f(x) = 0 (23) 
The solution of Eqs. (22) and (23) is 
f(@) = wo S [nan /(K + n%)] cos mt _ (24) 
n=1 


Eqs. (13), (17), (19), (21), and (24) now furnish the 
following evaluation of the integral Q: 


1946 


Q = (want/4) O K(aa™)*/(K + nt) (25) 


Hence, Eq. (16) is reduced to 
(a,™)? 


aida ene ye + Aim + 


| ( 4K Ai Mas bn 
l-vyp 


This is an explicit formula for the fixity that is developed 
by a wing structure that contains no supporting ele- 
ments for the compression panels other than longitudi- 
nal spars and rigid-acting bulkheads. 





(26) 





APPLICATION OF THE THEORY 


A number of numerical methods are available for 
evaluating the constant 2 and the Fourier coefficients 
a," which are required for the utilization of Eq. (26). 
An elementary procedure that is usually expeditious is 
illustrated below by an example. 

Values of ordinates z that locate centroids of elemen- 
tal plate-stringer areas with respect to the principal 
axis of inertia of the wing cross section are usually given 
in formal wing stress analyses. This tabulation defines, 
at discrete points, a functional relation between the di- 
mensionless coordinates § = rx/B and ¢ = 2/B. A 
fourth-degree polynomial usually represents this funce- 
tion with sufficient accuracy. Thus, from given wing 
data, the following relation has been derived by forming 
a table of divided differences and by applying Newton’s 
interpolation formula :? 
¢ = —0.00111£4 + 0.0172&* — 0.0970&? + 

0.173 + 0.0685 (27) 


The same relation could have been derived directly— 
though with greater labor—by assuming that a fourth- 
degree approximation , 


te 


€ = co&t + cr&® + co€? + sé + 0.0685 


is admissible, and then solving four linear equations that 
result by selecting four representative points from the 


given geometric data. 
Differentiation of Eq. (27) yields the following re- 


sults: 
af’ = —0.0139E* + 0.162? — 0.610 + 0.545 \ (28) 
ro" = —0.0417& + 0.3245€ — 0.610 
With the following formula of geometry, these equa- 
tions determine the curvature of the airfoil profile: 
= B/rR = —xf"/(1 + wg’)? (29) 
Since the effect of the denominator in this expression is 


not large, it may be expected that « is roughly propor- 
tionalto¢”. A plot of Eq. (29) reveals that a reasonable 


approximation is 


xk = —0.85x¢" = 0.0355¢? — 0.276 + 0.520 (30) 


The 
typical 
corresp 


= 1 


With tt 
the val 


correspo 
0.033, 0. 
Likewise 
are 1.10, 
denoting 
ling by ¢ 
1.64, C3 





(25) 


(26) 


oped 
 ele- 
tudi- 


e for 
ients 
(26). 
US is 


men- 
cipal 
viven 
fines, 
1e di- 
. OA 
funce- 
wing 
ming 
ton’s 


s that 
m the 


1g Te 


sion is 
ropor- 
ynable 


(30) 








STABILITY OF SEMIMONOCOQUE WING STRUCTURES 


Now, by the theory of Fourier series, the coefficients 
a," of Eq. (21) are represented by the integrals 


an) = (2/n) for sin mg sin ng d= (31) 


Hence, by Eq. (30), 


an!” = 0.0226 fo*é sin mé sin né dé — ; 
0. 70 f" & sin mé sin né dé 3 aie oe (32) 


m~AnN 
Also, by virtue of Eqs. (14) and (27), 


f = 2/B = —0.00071 fo"é sin? mé dé + 
0.0109 foé* sin? mé dé — 0.0617 fo" sin? mé dé + 
0.110, f,"& sin? mé dé + 0.0685 (33) 


Since the integrals of Eqs. (32) and (33) are generally 
useful, they are evaluated below: 
Je" sin mé sin né dé = | 
(—1)™ + mr { (m — nn)? — (m+ n)~*},m ~n 
So7é sin mé sin né dé = 1/2{1 — (—1)™ + n} Xx 
{(m +n)? — (m—n)-},m#n 


m 5 3 
Je sint me dg =F — Ph... 








10 2m? 4m4 q (34) 
4 2 
, & sin? mé dé = _ _ 3n° 
8 8m? 
r? 
sin? mé d = Sool 
ll £ si 1 dé = ; on 
i? & sin? mé di = 17/4 
With these results, Eqs. (32) and (33) yield 
a,” = 0.0710 [(m — n)-* — (m + n)-*], 
m #n,(—-1)™ +2 = +1 
a,” = 0.105 [(m — n)-? — (m + n)-*], (35) 
m #n,(—1)™ +" = —1 
a,” = 0.202 — breif /m?) 
e = mn 9 
-_ = + = 0800+ 0.116 0.0102 (36) 
Zmar. a8 m? m4 


The foregoing results were obtained from data for a 
typical wing, for which the pertinent constants (Eq. 10) 
corresponding to a station near the root are 


A = 10.83, A; = 1.17 X 10~-, A, = 1.28 X 10-5, 
K = 454 


With these constants, it is determined by Eq. (35) that 
the values of the expression 


E(ex™)*/[1 + (n*/K)] 
corresponding respectively to m = 1, 2, 3, 4, 5, 6, 7 are 
0.033, 0.0315, 0.0213, 0.0144, 0.0099, 0.0072, and 0.0053. 
Likewise, by Eq. (36), the respective values of Zmaz./Z 
are 1.10, 1.22, 1.25, 1.25, 1.25, 1.25, and 1.25. Hence, 
denoting the value of C corresponding to m-lobe buck- 
ling by C,,, one obtains, by Eq. (26), C, = 1.49, G = 
1.64, C; = 1.54, Cy = 1.46, Cs = 1.41, Cy = 1.38, and 
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Fic. 2. Typical cross section of test beam. 


C; = 1.38. Beyond C;, the numbers of this sequence 
obviously increase. Thus, the design value for the 
fixity is C = 1.38 in this example. 


TEsT RESULTS 


Five beams with cross sections like that shown in Fig. 
2 have been tested. The stringers were '/, by 5/s in. 
24S-T84 unclad bars spaced at 1!/, in., and the skins 
were 0.05l-gage or 0.080-gage 24S-T86 Alclad sheets. 
The bulkhead chord members were '/s by 11/4 by 1'/¢ 
in. plain 24S-T Alclad angles, with two intermediate 
supports as shown in Fig. 2. All specimens had con- 
stant radii of curvature on their upper surfaces. 

The specimens were tested by applying bending 
couples at the ends. Injurious end effects were pre- 
vented by local reinforcements extending over two bays. 
In each case, the unreinforced test section extended over 
five bays. The loading mechanism was operated by 
means of a calibrated system of hydraulic rams with 
pressure-gage attachments. The measured loads deter- 
mined the bending moments, and, since pure flexure was 





Fic. 3. 


Buckling of semimonocoque wing structures. 
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TABLE 1 
Descriptive Constants and Calculated Characteristics of Test Beams 

Beam No. 1 2 3 4 5 
L 8 8 8 12 12 
B 30 30.1 31.4 31.4 31.4 
R 2 73.5 30 30 30 
t 0.050 0.051 0.051 0.051 0.080 
be 0.173 0.172 0.173 0.167 0.196 
p 0.089 0.089 0.089 0.0915 0.0974 
kK 0 0.131 0.333 0.333 0.333 
K 4.92 4.96 5.40 2.40 2.40 
A 18.5 18.6 17:3 85.5 101.0 
A, 4.20 X 10-5 4.40 X 1075 3.70 X 1075 1.85 X 1074 5.39 XK 1074 
Az 0.00730 0.00717 0.00665 0.00976 0.00735 
m 1 3 5 5 5 
Zmaz. /3 1.00 1.05 1.14 1.13 1.13 
c 1.01 1.25 1.78 2.57 3.04 
en 13,000 16,100 22,800 15,500 20,800 





applied, these bore a known ratio to the plate-stringer 
In some cases the stresses were checked with 
On all beams, deflec- 
Fig. 3 


stresses. 
electric strain gage equipment. 
tions were thoroughly mapped with dial gages. 
is a photograph of a typical test setup. 

Dimensions of the specimens are shown in Table 1. 
Aside from differences indicated by these dimensions, 
the specimens were identical. The following material 
constants were adopted for all specimens: 


E = 10.5 X 108 Ibs. per sq.in. 
v = 0.30, A = 1.00 


The torsional rigidity J of a bulkhead chord member 
is obtained from the well-known torsion formula for 
open sections: 3J = Apt”. For the '/s by 11/4 by 
1'/, in. angle sections used in the box beams, the area 
A, is 0.280 sq.in., and, consequently, J = 0.00146. 

The constants A, Ai, As, K, and x given in Table 1 
were computed by the formulas of Eq. (10). The ordi- 
nates 2 were measured from the centroidal axes of the 
beam cross sections, and the quantities 2 were obtained 
by numerical integration of Eq. (14). Since all beams 
had circular camber, the Fourier analysis (Eq. (21)) de- 
generates, in each case, to 


a,” = kiim=n 
” Oifm ~n 


Hence, Eq. (26) becomes 





Z Ax? 
Ce Fe EB mss hg 
(=) ‘Tree 
(=4 4 Az)m* (37) 
l-—yp 


From this formula, the integer m that furnishes the 
smallest value to C has, in each case, been determined by 
trial. Values of m and C so obtained are given in Table 
1. Finally, the anticipated buckling stresses have been 
computed by the Euler column formula o,, = Cr*E X 


(p/L)?. 


It may be observed that the effects of the bulkhead 
uprights have been neglected in the calculations of 
Table 1. For Beam | these effects are certainly negli- 
gible, and for Beam 2 they are null because the nodes 
in the twisted bulkhead chord occur at the uprights. 
For Beams 3, 4, and 5, the effects can be no greater 
than those attending an increase of m from 5 to 6, since, 
with m = 6, there is no bending of the uprights. Hence, 
for the first four beams, the perturbation of C due to 
the uprights is not more than 1 per cent, while for Beam 
5 it is known only that it is less than 9 per cent. 

To substantiate the foregoing theory, it would be 
simplest to compare the calculated buckling stresses 
(Table 1) with observed buckling stresses. However, 
because of unavoidable imperfections of materials and 
workmanship and because of natural deformations at- 
tending bending of the beams, ideal buckling did not 
occur. A picture of the true behavior of the specimens 
is given by stress-deflection curves (Fig. 4). These 
show deflections relative to the spars at or near the 
centers of the compression surfaces. Relative deflec- 
tions were automatically determined by mounting the 
dial gages on chordwise brackets that rested on the 
spar flanges (Fig. 3). Thus the curves (Fig. 4) present 
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a true representation of the chordwise arching or hump- 
ing of the compression panels. 

The deflection curves (Fig. 4) show definite knees, 
which presumably lie near the ideal buckling stresses. 
The calculated buckling stresses (Table 1) are indicated 
on the curves by stars, and the corresponding fixity 
factors are appended. It is seen that, in each case, the 
calculated buckling stress lies on the knee of the corre- 
sponding stress-deflection curve. This circumstance is 
a good confirmation of the present theory. 


CONCLUSIONS 


It is apparent from the test results illustrated in Fig. 
4 that the conventional practice of utilizing an invari- 
able specified fixity factor for all wing stress analyses 
admits significant refinements. The frequent specifi- 
cation C = 1.5is certainly unsafe if the wing surface is 
nearly flat in the interspar region, whereas it is exceed- 
ingly conservative under other conditions. 

Since the factor C determines reduced column lengths, 
its influence is smaller in the short column range than in 
the Euler range. However, with the advent of pre- 
cipitation heat-treated aluminum alloys, column design 
curves have become much steeper, and fixity in the 
short-column range has consequently assumed greater 
importance than formerly. 

The foregoing theory is presented as a rational ap- 
proach to the “fixity problem.”’ Within the range of 
the experiments, the theory has been well confirmed. 
Further tests are needed to verify the theory in cases of 
noncircular camber and to determine the effect of pre- 
mature buckling of the skin. Although all of the test 


specimens carried appreciable loads in excess of buck- 
ling, it is probable that this effect has little significance 
when the stringers are thin-walled sections that cripple 
under excessive bending. 

The fact that the fixity due to membraiue action of 
the skin varies as the fourth power of the bulkhead 
spacing (cf. Eqs. 10 and 26) presents the possibility of 
designing wings with comparatively few bulkheads. 
However, in order to realize high fixity from membrane 
action, it is necessary to limit the buckling pattern to 
modes that exhibit a small number of chordwise waves— 
preferably to the fundamental mode defined by m = 1. 
This is most directly accomplished by introducing, 
between the bulkheads, chordwise ribs that exert small 
restraints against single-lobe buckling but which never- 
theless control the buckling pattern. Results obtained 
by generalizing the foregoing theory indicate that ex- 
ceedingly high fixities may sometimes be obtained with 
small weight cost by this device. However, test results 
for this type of construction are lacking at present. It 
should be observed that the theory, as presented, is not 
applicable for cases in which C exceeds four, since only 
single-loop buckling in the spanwise direction has been 
contemplated. 
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Reduction of the Adhesion of Ice to 
De-Icer Surfaces” 


D. L. LOUGHBOROUGH}? anv E. G. HAASt 
The B. F. Goodrich Company 


SUMMARY 


Ice is removed from a De-Icer surface by mechanically applying 
concentrated stresses to the ice-rubber interface. It is obvious 
that a low value for the adhesive bond should be desirable, but it is 
conceivable that a zero value might not be practical. The shear 
strength of the bond to the De-Icer surface is as high as that be- 
tween ice and any other solid; it increases linearly with decrease 
in temperature from a value of zero at 0°C. to 150 Ibs. per sq.in. 
or more at —25°C. It is approximately 15 per cent higher to the 
surface of a rubber that has soaked in water and probably de- 
pends on surface smoothness only to the extent that a rough sur- 
face presents a greater area. 

It has been found that silicone, either formed in situ by poly- 
merization on the De-Icer surface or applied as ready-made poly- 
mer, improves the performance of De-Icers by lowering the ad- 
hesion of ice to the surface by as much as 90 per cent. These 
treatments have shown good results both in the laboratory and in 
flight tests. It is expected that, because of the toughness of the 
films formed by silicones, they will need to be applied only a 
few times during an icing season. 


INTRODUCTION—PRESENT METHOD OF ICE REMOVAL 


bia METHOD for the removal of ice from aircraft 
that has been used most extensively in this country 
is that involving the use of a pneumatically operated 
rubber boot known commercially as a De-Icer. The 
De-fcer consists of a rubber sheet, containing cells or 
tubes, stretched and mechanically attached at its outer 
margins to the leading-edge surface to be protected. 
The inflation and deflation of the cells cracks the ice so 
that it is carried away by the air stream. Fig. 1 isa 
schematic drawing of the basic three-tube De-Icer in 
the three inflation stages, while Fig. 2 shows an actual 
installation. 

Though this device works well under most icing condi- 
tions, difficulty sometimes is encountered in completely 
removing ice, especially at low temperature where the 
ice adheres particularly tenaciously. For some time it 
has been felt that performance would be improved if 
the adhesion of ice to the surface could be lowered. 
The work toward this end reported ‘here is only one 
phase of the research and development work that is 
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being followed to improve the performance of this prod- 
uct. 

Although the De-Icer is constructed of a pure gum 
natural rubber stock with fabric reinforcement where 
necessary, the surface usually is a neoprene cement 
coat. When the De-Icer surface is mentioned it is this 
thin neoprene surface which is meant. 





STRETCH AREA 
TUBE AREA 


CENTER TUBE 
INFLATED 






INFLATED 


_ fn OUTER TUBES 


Schematic drawing of basic three-tube De-Icer in three 
inflation stages. 


Fic. 1. 





Fic. 2. De-Icer installed on horizontal tail of an airplane. 
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REDUCTION 


PHYSICAL PRINCIPLES UNDERLYING ADHESION 


A full understanding of the role that the adhesion of 
ice to rubber plays in De-Icer operation requires a con- 
sideration of some of the physical principles involved 
in the operation. A detailed discussion of these matters 
is complicated and outside the scope of this paper, but 
certain fundamental matters will be considered. 


Applied Stress 


On applying suitable stresses to a rubber sheet cov- 
ered with ice the two components separate. The type 
of stress applied is not critical and is determined solely 
by mechanical convenience. In the De-Icer, forces in 
two different directions are used: (1) in the tube area 
(Fig. 1) the interface is distorted in a direction perpen- 
dicular to its plane; (2) in the stretch area (Fig. 1) the 
interface is distorted in a direction parallel to its plane. 

The applied stress is magnified by the stress concen- 
tration to such an extent that the local force on the in- 
terface exceeds the bond strength. Analyses of these 
forces have been made and show that the multiplication 
of the stress is many fold. The bond can then be 
broken with a low load applied slowly. This is the im- 
portant fact about De-Icer operation. A thin steel 
shell could be made to remove ice in the same way the 
rubber boot does, except extremely high forces would be 
required. The interface distorts until the concentrated 
force exceeds the adhesive force, and then the ice film 
cracks loose from the surface. The object of De-Icer 
design is to make the necessary applied load as low as 
possible. This can be done, as explained above, (a) by 
using a low modulus, elastic material so that high stress 
concentrations are possible, and (b) by lowering the 
adhesive force. It should be pointed out, however, 
that the latter scheme must not be carried too far, for, 
if the adhesive force were lowered to a zero value, it is 
doubtful that the ice cap could be removed. The stress 
applied by the De-Icer could not be transmitted to the 
ice to cause the breaking of the cap. It is desired, then, 
to reduce the adhesion to some low nonzero value. 


Stress Concentration 


The measurement of the adhesive force between two 
solids is complicated by the fact that stress on the inter- 
face is concentrated and locally exceeds the applied 
stress. When the materials are as different as ice and 
tubber, all values of adhesion must be questioned, 
since results from two slightly different experimental 
atrangements may easily differ by a factor of five. It is 
important for the evaluation of the true adhesion that 
the applied stress be constant over the interface or that 
the stress distribution be accurately known. 


Methods of Measuring Adhesion 


Each of the many methods which have been used to 
measure the adhesion between ice and various solids has 
an advantage for particular conditions. A centrifugal 
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method developed by the Massachusetts Institute of 
Technology De-Icing Research Laboratory has been 
used to a large extent in this work. The apparatus 
(Fig. 3) consists of a horizontal plate that can be ro- 
tated at controlled speed. An ice button is frozen to a 
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Fic. 3. Apparatus for measurement of ice adhesion by,centrif- 
ugal method. 


Rotating plate. 
High-speed motor. 
Sample. 

Ice button. 

Guard ring. 
Window. 


Section AA—Enlarged cross section of sample and ice button 
during measurement to adhesion. 
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Fic. 4. Effect of thickness of sample on apparent adhesion of 
ice to pure gum rubber, —25°C. 
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rubber sample cemented to the disc, and the centrifu- 
gal force required to just free the button is measured. 
The motor B drives the plate A at the desired speed, 
measured by a stroboscope through the window F of the 
refrigerator. The speed is increased until the ice button 
D flies off the test sample C and strikes the guard ring 
E. In Fig. 3, Section AA, is shown an exaggerated de- 
tail of the ice button just prior to removal. Since the 
ice is: above the plane of rotation, a couple is developed 
which rotates the sample about an axis parallel to the 
plane of rotation and at right angles to the radius of the 
plate. The applied centrifugal force is not, then, a pure 
shear unless the test piece (C) becomes very thin or 
very, hard, Fig. 4 gives the values measured at various 
thicknesses and shows clearly that the apparent ad- 
hesion is a function of the undistorted rubber sample 
thickness as indicated by Fig. 3, Section AA. 

Another device for measuring the adhesion is shown 
in Fig. 5. In this experimental arrangement the inter- 
face is subjected to almost pure shear by rotating brass 
shells A and B with respect to each other. Values se- 
cured by this technique are given in Table 1. 





TABLE 1 
Torsional Method of Measuring Adhesion 
Mechanical Goods Rubber Stock 





Temperature Adhesion 
(“<) Lbs. per Sq.In. 
—10 80 
— 20 170 
—30 200 





Other methods of determining the adhesion based on 
totally different principles have been used. Two are 
given below. 

(1) Calculation of the stress concentration around a 
piece of ice adhered to a large piece of rubber in tension 
(as in the stretch area) has been made. The necessary 
mathematics for this calculation was developed by 


























Fic. 5. Apparatus for torsional method of measuring adhesion. 


A. External brass shell. 
B. Internal brass shell. 
C. Rubber. 

D. Ice. 
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Donnell.1_ By measuring the equilibrium area of con- 
tact between the two solids when a fixed elongation is 
applied to the rubber, it is possible to get an estimate of 
the adhesion. The experiment is too difficult to con- 
trol for extensive use and has value primarily because it 
indicates the mechanism of De-Icer operation. No de- 
tails will be given. 

(2) By calculating the shear stresses between rub- 
ber and ice when the interface is deformed perpendicular 
to the interface, it is possible to secure an approximate 
value of the adhesive force. 

In all of these methods the ice formed is of the hard 
glaze type. No method of easily forming other types is 
available in our laboratory at the present time. It is 
probable that rime or any other air-containing ice will 
show a lower value of the apparent adhesion only be- 
cause of the reduced contact area. 

Because the interfacial area is not accurately known, 
only a few results are given in absolute units. In most 
instances the values listed are comparative, and indi- 
cate either a per cent of the value for the untreated sur- 
face or a per cent of the value for a pure gum stock. 


True Value of the Adhesion 


From the first two tests it can be seen (Table 1, Fig. 
4) that the adhesion of ice to rubber-like materials is 
about 170 Ibs. per sq.in. at —25°C. In the second 
method discussed the difficulty of preparing samples is 
too great for routine testing. In general, the centrifu- 
gal device has been chosen, using test pieces about 
0.05 to 0.1 in. thick. Although it was realized that the 
absolute adhesion was not correct, it was believed that 
the relative values for the different surfaces would be. 
Where some notion of the absolute value of the ad- 
hesion of ice to flexible material was desired, thin films 
were used as test samples. 

That the true adhesion of ice to rubber is as high as 
that of ice to many other materials is shown in 
Table 2. 





TABLE 2 
Adhesion of Ice to Miscellaneous Materials 
—25°C.—Centrifugal Method 


Shear Strength of 
Bond (Lbs. per Sq.In.) 





Test Material 


Ice 250 
Steel—polished 185 
Copper—polished 124 
24ST aluminum 220 
Glass (96 per cent silica) 158 
Microscope slide 122 
Polyethylene 160 
Rubber 170 





From the above discussioti it follows that the use of 
rubber in De-Icers is not dictated by the low adhesion 
but by the fact that low applied stresses can be used 
because of the stress concentration. The high value of 
the adhesion does indicate the. possibility of consider 
able reduction, but, before discussing’ methods of 
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achieving this, some understanding of the factors that 
affect the magnitude of the adhesion should be help- 
ful. 


Factors Affecting the Adhesion 


Temperature-—The adhesion of ice to rubber and 
other materials changes linearly with the temperature 
as shown in Fig. 6. The slope of the curve for the thin 
sample is over four times as great as that for the thick 
sample. The distortion resulting from the 0.1-in. 
thick rubber sheet reduces the adhesion 76 per cent at 
all temperatures. The values for the thin sample ap- 
proach more closely the true bond strength between ice 
and rubber. 

The strong dependence of adhesion on temperature 
is not restricted to that between rubber and ice. As is 
shown in Table 3, the slopes are high for all materials 
but two. One of these is for the ‘‘adhesion”’ of ice to ice 
and forms an interesting and unexplained exception to 
the other materials. The other exception is silicone 
on rubber, which will be discussed later. 





TABLE 3 
Temperature Dependence of Adhesion 





Temperature Slope 


Test Material (Lbs. per Sq.In. per °C.) 
0 


Ice 

Steel —7.4 
Copper —5.0 
Aluminum —8.8 
Glass (96 per cent silica) —4.5 
Polyethylene —-5.8 
Rubber —6.0 
Silicone treated rubber —0.4 





This striking dependence on temperature (6 lbs. per 
sq.in. per °C.) is of interest to those studying the theory 
of adhesion. It is of even more importance for the 
understanding of the mechanical removal of ice at low 
temperature. It is under these conditions that a re- 
duction in adhesion can be expected to give the greatest 
improvement in performance. 

Surface Roughness.—The lack of sufficiently detailed 
information on molecular forces has prevented the de- 
velopment of any theory regarding the origin of the 
adhesion between two solids. It has generally been as- 
sumed that some materials adhere by the interlocking of 
surface irregularities, though the proof of this has never 
been secure because of lack of knowledge of the true 
contact area. In general, most workers imply that the 
adhesion largely results from molecular interaction. 
This is particularly probable when the adhesion is 
strongly temperature dependent in a temperature range 
in which neither solid greatly alters physical proper- 
ties. 

Two sets of data are presented in Table 4 to show the 
dependence of ice adhesion on surface conditions. In 
Part A of the table the rubber was cured against sur- 
faces of different visually estimated roughnesses. In 
Part B the rubber was cured against a smooth surface 
and scored in the ways described. 


Surface Roughness 


—30°C.—Pure Gum Rubber 





A Apparent 

Adhesion 

Visual Per Cent 

Surface Cured Smoothness of Highest 

Against Rating Value 
Fabric—rough 6 100 
Holland cloth backed 

by fabric 5 50 
Parchment paper 4 70 
Holland cloth 3 60 
Aluminum 2 65 
Glass 1 25 





— 20°C.—Neoprene Surface 
B Adhesion (Per Cent 
of Value for the 


Description of Sample Original Surface) 


Pricked with a pin 125 
Covered with Santocel 164 
Cut with razor blade 182 
Roughened with emery 210 





The increase in adhesion resulting from roughening 
of the surface could probably be completely attributed 
to the increase in contact area. This has not been 
proved. No matter what explanation is finally given, 
it is obvious that it is important to maintain the 
smoothness of De-Icer surfaces. 

Moisture Content——The data in Table 5 show the 
ratio of the adhesion value of dry to wet material for 
several different .materials. In making the measure- 
ments, the material that had recently been soaked in 
water was considered to have a wet surface even though 
the surface was wiped dry. 
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TABLE 5 
Effect of Moisture on Ice Adhesion 
Soaked in Water at 70°C. for 1 Hour 





Ratio of Adhesion Value 


Rubber Stock of Dry to Wet Material 


Pure gum rubber 0.83 
Rubber plus 60 parts Gastex 0.91 
Neoprene compound 0.92 
De-Icer surface 0.62 
Mean of measurement on 10 

different stocks 0.87 





In De-Icer operation the adhesion is almost always 
to a wet surface. The water absorption of rubber is 
not primarily due to the rubber hydrocarbon but to 
the pigments, protein impurities, etc., which are present. 
Samples of gum rubber prepared from latices specially 
washed to remove moisture-absorbing ingredients were 
found to have only about half of the adhesion of a nor- 
mal pure gum rubber stock. 

Tables 1 through 5 have indicated the nature of ice 
adhesion. The desirability of finding a way to reduce 
it is obvious. 

Surface Treatments——The adhesion of ice to most 
elastomers from which De-Icers can be made is so 
nearly the same that little advantage could result from 
changing the elastomer. The only known exception to 
this generality is silicone ‘‘rubber,” which shows prom- 
ise especially for use as a surface coating on a rubber 
stock. 

Present De-Icers have a coating of neoprene which 
protects the rubber from the sunlight and, because of 
its compounding ingredients, renders the surface con- 
ducting. Thus it is protected from electrical dis- 
charges. This surface cannot easily be altered unless 
the new material can also fulfill these functions. As an 
example Table 6 presents the adhesion of ice to a num- 


ber of elastomers. 





TABLE 6 
Adhesion to Thin Films of Various Elastomers 
—25°C., 0.002 In. Thick—Centrifugal Method 
Adhesion 
(Lbs. per Sq.In.) 





Elastomer 
Rubber 170 
Hycar OR 150-160 
GR-S 130 
Thiokol FA 138 
Koroseal! 120 
Neoprene? 150 
Polyethylene 160 
Polystyrene® i15 
Silicone rubber‘ 45 





1 Value determined by extrapolation. The value measured, 
—25°C., is certainly too low, because the film is weak and ad- 
hesion of the film to metal (surface) is poor. 

2 For comparison purposes this may be considered as a De-Icer 
surface. ’ 

3 Not an elastomer. Included here for general interest. 

4 Estimated from a determination made on a thick piece. 
Value on a 0.1-in. sample was 10 Ibs. per sq.in. (See Fig. 6.) 


A great number of surfaces and surface treatments 
have been tried by this and other laboratories in the 
past with little success. In a recent incomplete list of 


materials which had been tried, over 1,000 distinctly 
different chemicals were recorded. 

In general, those materials that have been found to 
increase the ease of ice removal have done’so without 
decreasing the ice adhesion. Three general types may 
be listed: (1) a weak salve or grease that breaks and 
comes off with the ice; (2) a material that is a freezing 
point depressant; and (3) materials that decrease the 
contact area. When these mechanisms are used in 
practice, they have a short life or require continuous 
replenishing during time of greatest need. 

Despite the long record of failures, it appeared that 
the most promising approach to the problem of De-Icer 
improvement was through the use of surface treatments, 
There are some who claim that a low-adhesion surface 
must be hydrophilic, perhaps as a result of confusing 
the effect of the hydrophilic nature of a treatment with 
the freezing point depressant effect. All of the low- 
adhesion surfaces that we have tested have been hydro- 
phobic, but not all of the hydrophobic surfaces have 
exhibited low adhesion. 


SURFACE TREATMENT WITH SILANES AND SILICONES* 


Of the number of materials tried, silicone, either 
polymerized from silanes applied to the test surface or 
applied directly as a silicone, proved most promising 
for use on De-Icers. Detailed investigation shows that 
they almost completely fulfill the requirements. The 
balance of this paper will be devoted to a discussion of 
the results with these materials. 


Silanes 


The silanes referred to in this work are organo-silicon 
halides analogous to the substituted methanes and 


having a general formula 
R,Si Clin 


On exposure to a moist atmosphere they react rapidly 
with water-forming silicols. The silicols polymerize to 
silicones at a rate depending on the nature of the hydro- 
carbon radical. The dimethyl derivative polymerizes 


rapidly to an oil. 
(CH3)2SiCl, + HO — (CHs3)2Si(OH)2 + 2HC1 
CH; CH; CH; 


— —O—Si—O—Si—_O—Si—0 + H;0 


| | | 


CH; CH; CH; 


The chemistry and properties of these materials have 


-been described in the literature.?~* 


. Materials of this nature and for this purpose are supplied 
commercially by The B. F. Goodrich Company. Their use for 
lowering of ice adhesion is now the subject matter of patent 


applications. 
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In Table 7 are given the results of some tests on the 
adhesion lowering resulting from treating rubber sur- 
faces with liquid silanes, which, of course, polymerize to 
silicones. 





TABLE 7 
Adhesion Reduction by Silanes 
— 25°C.—De-Icer Surface 





Adhesion 

Silane (Per Cent of Un- 
(Applied as the Pure Liquid) treated Surface) 
Monomethy] trichlorsilane 12 
Dimethyl] dichlorsilane 20 
Trimethyl monochlorsilane 23 
Mostly diethyl dichlorsilane 2 
Mostly dibutyl dichlorsilane 25 
Mostly diphenyl] dichlorsilane 25 
Phenyl butyl dichlorsilane 17 





In this type of testing it is difficult to control the ex- 
periment completely. Here and in subsequent tables, 
differences in adhesion values of 5 per cent are not sig- 
nificant. 

Other methods of applying the material are also satis- 
factory. 








TABLE 8 
Method of Applying Silane 
— 25°C.—Silane 
Adhesion 
(Per Cent of 
Surface to Untreated 
Treatment Which Applied Surface) 
Vapor in a closed container : 
50 min. Rubber 51 
Vapor in a closed container 
100 min. Meera 30 


Soaked in water 1 hour at 
70°C., then exposed to 
vapor 90 min. ae 40 

Equilibrium with room hu- 
midity, then exposed to 








vapor 90 min. sisi 45 
Hung over liquid—10 min. De-Icer surface 18 
Hung over liquid—60 min. De-Icer surface 15 
Treated with liquid De-Icer surface 12 
Treated with solutions See Table 9 

TABLE 9 


Applying Silane from a Solution 
A. Various Concentrations in CCl, on De-Icer Surface 





Concentrations Adhesion (Per Cent of 
(Per Cent by Volume) Untreated Surface) 

4.5 . 

9 33 

17 22 

28 16 

33 23 

50 18 

100 12 





B. Various Solvents 


30 Per Cent 50 Per Cent 

Solvent by Volume by Volume 
Naphtha B 37 ~ 
CCl 23 18 
CeHuy 42 21 
C.H,Ch 33 15 
CeHe 33 9 
CsH;Cl 49 36 





The adhesion is lowered by all methods of treatment, 
but the use of the pure liquid gives slightly better re- 
sults. In most of our tests this method has been used 
exclusively. 

A successful film must, in addition to lowering the ad- 
hesion, be reasonably permanent. Though the durabil- 
ity of the film will be considered in detail later, the data 
that follow will give the adhesion results after expo- 
sure to an ultraviolet light and then to a weatherom- 
eter. 

The ultraviolet light is a Hanovia Chemical Company 
lamp, Type 7420, placed 15 in. from the sample in a 
partially closed space. The temperature rose about 
15°C. above room temperature. The weatherometer is 
similar to that used in testing painted surfaces. The 
samples are suspended in a circle about a carbon arc 
light and rotated about the light every 15 min. Once 
during each rotation they are subjected to cold-water 
spray. 

In Table 10 are collected the results achieved by using 
the pure silanes. The differences among these treat- 
ments show most significantly after aging under an 
ultraviolet light. 





TABLE 10 
Adhesion Reduction by Certain Silanes 
Pure Liquids 





-———Adhesion (Per Cent of Control) —— 
After 86 Hours 
Under Ultraviolet 


After 86 Hours Light and 20 





Under Ultra- Hours in 
Silane Initial violet Light Weatherometer 
CH;SiCl; 12 >100 % 
(CHs)2SiCl, 16 41 62 
(CHs)3SiCl 23 75 a 
Silicones 


Two possible disadvantages in polymerizing the sili- 
cones in situ are apparent: (1) the polymerization con- 
ditions are uncontrolled, and (2) the copious evolution 
of HCl may be hazardous to personnel and equipment 
unless special precautions are taken. It would be more 
convenient to make the silicones and apply them to the 
surface to be protected, provided the results were 
equally good. Which results are better may depend 
on how firmly the respective films are bound to the rub- 
ber surface. It is possible that the mode of attach- 
ment of the rubber to the film formed from the silanes 
is through a primary valence brought about by the reac- 
tion of the halogen in the silane with a hydrogen in the 
surface. Table 11 and Fig. 7 present the data on sili- 
cones and mixtures of silicones. The adhesion reduc- 
tion is noticeably improved by lowering the CH;-Si 
ratio to a value of 2 or less. The results do not depend 
on the silanes chosen but on the methyl silicon ratio. 
Similar results have been found for mixtures of silanes 
except that in this case copolymers of mono- and tri- 
methyl silane are significantly poorer than other copoly- 
mers of the same CH;/Si. 
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e MONOMETHYL TRICHLORSILANE & 
DIMETHYL DICHLOR SILANE. 

o DIMETHYL DICHLORSILANE & 
TRIMETHYL MONOCHLOR SILANE. 

eMONOMETHYL TRICHLORSILANE€. 
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Fic. 7. Adhesion to De-Icer surface treated with silicones 


prepared from mixtures of methyl chlorsilanes, —25°C. Thick- 
ness of sample, 0.060 in. 














TABLE 11 
Silicones from Mixtures of Silanes 
Adhesion 
—-(Per Cent Untreated Surface) 
86 Hours in 
Ultraviolet 
86 Hours’ Light and 20 
Under Hours in 
Ultraviolet Weatherom- 
Monomers CH;/Si Unaged Light eter 
Mono- and di- 
methyl 1.75 6 40 50 
1.88 14 42 52 
1.90 13 62 54 
Di- and tri- 
methyl 2.08 8 52 65 
2.10 13 66 77 
2.20 16 73 te 
2.50 48 81 79 
Mono- and tri- 
methyl 1.67 17 48 7. 
2.00 37 85 87 
2.30 . 40 77 67 
2.50 42 81 81 
2.80 31 87 71 





The silicones can be applied in a number of ways, but 
a direct application is the most successful. 

The effectiveness of the silicones for the purpose is 
remarkable, especially since other materials, such as 
polyethylene, which present an equally pure hydrocar- 
bon surface, are poor for lowering ice adhesion. Further 
increasing the length of the aliphatic chain in the mono- 
meric silane does not improve the ice-adhesion lowering. 
Some other feature of molecular structure must account 
for the unique value of this polymer. The explanation 
finally accepted must account for the low ice adhesion 
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Fic. 8. Laboratory test of ice removal from treated and un- 
treated De-Icer surface. Three-tube De-Icer, —25°C. Ice 
thickness, 0.030 in. Methyl chlorsilane treatment. Treated 
surface on right. 


in a way that does not require high water repellency, 
for it can easily be shown that there is poor correlation 
between these two properties. 


Durability 


That lowering the adhesion improves De-Icer 
performance is indicated in Fig. 8. A laboratory 
wing section was mounted in a refrigerator with the 
leading edge up. The area on the right was treated with 
Icex prior to coating the rubber with a 0.030-in. thick 
cover of glaze ice at —32°C. The De-Icer was inflated 
once, and all of the ice came off the treated area. In 
the photograph the black areas are free of ice; in the 
white areas the ice is free but is held in piace by that 
which is tightly adhered and appears dark gray. This 
type of testing indicates promising materials but fails 
as a service test because it does not subject the surface 
to the severe abrasion met in service. To duplicate such 
action in the laboratory is difficult. 

The films, when placed on De-Icers in service, must 
withstand ultraviolet light, exposure to impact of water 
drops, and repeated removal of ice caps. The duplica- 
tion of these conditions in the laboratory is not strictly 
possible, but from certain of the tests described below 
some notion of the performance characteristics of the 
material can be gathered. Neither material has a meas- 
urable effect on the usually measured propertics of rub- 
ber or De-Icer surface, such as modulus, tensile strength, 
elongation, flex life, etc. 
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TABLE 12 
Durability to Solvents 
Silane 
Per Cent of 
Surface with 
No Silane 
Washed with benzene 4 
Washed with acetone 19 
Washed with butyl acetate 44 
Washed with methyl] ethyl ketone 33 
Washed with soap and water 19 
Soaked in water 1 hour at 70°C. 15 
None 12 
a os 
Stability Under Ice-Forming Conditions 
Silane 


Repeat Values on Same Spot on Sample 
——Adhesion (Per Cent of Control)— 





Rubber Aluminum Steel 
lst button 17 19 16 
2nd button 12 28 98 
3rd button 15 74 vi 


4th button 21 


The treatment is surprisingly stable to a range of 
solvents and is not dissolved or freed from the surface 
by water. 

In De-Icer operation many ice caps are removed in 
succession. If the material is low in strength (like 
Vaseline), some of it will come off with each ice cap and 
consequently will appear satisfactory on one test but 
poorer on each successive test. The silane-treated rub- 
ber surface is satisfactory but a similar treatment on 
metal is not. q 

The preprepared silicones under some circumstances 
are more resistant to weathering conditions than sur- 
faces treated with silanes to produce silicones in situ. 
Silicones can be produced to have almost any viscosity 
between that of water and that of a heavy grease or 
rubbery material. These latter materials are of special 
interest from the standpoint of durability. On initial 
application they reduce the ‘‘adhesion”’ greatly—like 
any grease. Actually, the material is removed by the ice 
cap. With successive ice-button removal the adhesion 
rises to a constant value. At this point the film is thin 
and strong, and, henceforth, no further silicone is re- 
moved by the ice. Data for one such grease are given 
in Table 14. 

That this increase in adhesion was accompanied by 
removal of silicone was proved by adding a small por- 
tion of radioactive phosphorus to the silicone and fol- 
lowing each adhesion test with a measurement of the 
activity of the treatment with a Geiger counter. 

When a surface treated with a high-viscosity silicone 
is exposed to a high-velocity water spray (100 m.p.h.*) 


* The velocity cannot be directly compared to airplane speeds 
because the radius of the sample while being sprayed is smaller 
than the wing (2-in. diameter cylinder), and, hence, impact veloci- 
ties are higher in the laboratory. Further, the drops are larger 
and liquid water content of the spray is higher than in natural 
fogs. 
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TABLE 14 
Durability Under Ice-Forming Conditions 
High Viscosity Silicone Grease 
Repeated Values on the Same Spot on Sample 

















Number of Adhesion 
Removals (Per Cent of Control) 
1 3.6 
2 4.0 
3 5.0 
4 7.2 
5 7.8 
6 10.4 
7 14.0 
Ss 15.6 
9 17.0 
10 17.2 
11 13.4 
12 19.2 
13 19.8 
TABLE 15 
Durability of a Silicone Grease to Water Spray 
Adhesion 


—(Per Cent Untreated, Unexposed Surface). 


Hours Silane Silicone Liquid Silicone Grease 
0 
1 17 22 6.0 
2 80 - 
3 71 : 
4 87 85 
8 a ‘t 
16 vs = by 
24 $f “ 15.0 
48 28.0 
72 56.0 
96 81.0 











Effect of Stretching on Efficacy of Silane Treatment 


Adhesion (Per Cent of 
—Untreated Surface)— 





Treated Treated 
Exposure Stretched Unstretched 

No exposure 23 ‘ 19 
Stretched 30 per cent in weatherometer 

50 hours 58 55 
Stretched 30 per cent on roof—2 weeks 41 28 
Looped on roof—2 weeks 26 19 
Flexed (200 cycles per min.) on roof—2 

weeks 37 31 





the adhesion rises to the same value as indicated in 
Table 14 and remains constant during many hours of 
exposure. This is shown in Table 15 in comparison 
with a silane treatment. 

Since the De-Icer must stretch to operate the treat- 
ment must also elongate. Table 16 indicates that the 
silane-treated surface is satisfactory in this regard. 
The value after 2 weeks on the roof is the same whether 
the sample was flexed or not. 

Table 17 shows some typical weathering data on vari- 
ous films under a number of conditions. The water- 
spray data were taken as follows. High-velocity water 
spray was thrown on the sample in the refrigerator. 
About 0.1 in. of ice per min. was built on the surface. 
By the inflatable De-Icer principle the ice caps were 
removed approximately every minute. The adhesion 
values were measured at the end of the times noted. 








134 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 











TABLE 17 
Weathering Conditions on Various Films 
Adhesion (Lbs. per Sq.In.)————. 
Exposed on Exposed o 
Exposed to Roof in Roof in 
Material —Room, Wks.——California, Wks.——Akron, Wks.— 
vr a 6 a > as 
A silane hak ic A 25 47 #47 71 7. 2 
A silicone 244 7 8 24 26 37 14 22 33 





Exposed 4 Hours to High-Velocity Water Spray After Akron 
Roof Exposure of Following Times in Weeks 


Original 0 1 Wk. 2 Wks. 4 Wks. 
A silane’ 6 17 26 27 61 
A silicone P 40 34 40 32 





This table and those previously given indicate that 
the adhesion lowering accomplished by these com- 
pounds is reduced (1) by exposure to ultraviolet light 
at about 90°F., (2) by exposure to ultraviolet light and 
high-velocity water spray, and (3) by exposure to sum- 
mer weather. Further, it is not affected by exposure 
to normal room conditions. The serious damage done 
by the spray is nearly independent of the exposure to the 
ultraviolet light and has recently been shown to be 
equally detrimental when used at room temperature. 
Since the correlation between these laboratory condi- 
tions and those met in flight are not known, the labora- 
tory tests can only be used in a comparative way. 

In Fig. 9 is a photograph illustrating the ice-removal 
ability of the treatment after exposure to natural icing 
conditions as encountered on the summit of Mt. Wash- 
ington. This photograph is one of a great number se- 
cured in our extensive test program on De-Icer design 
and surface treatments carried out on the summit of 
Mt. Washington where the temperature, wind velocity, 
and humidity are such as to give frequent icing condi- 
tions comparable to those encountered by low-speed air- 
planes in flight. The tests were made with cooperation 
of the Weather Bureau staff assigned to the Mt. Wash- 
ington Observatory. The results of these tests prove 
that surfaces treated with silicone materials leave some 
50 per cent less residual ice. A silane was effective for 
over 15 days of almost continual ice and high winds. 
Silicone oil had an equal life. A specially treated sili- 
cone grease was still effective after 1 month of exposure. 

The results obtained in the laboratory have been con- 
firmed in flight tests by comparing the operation of 
treated and untreated De-Icers on airplanes in service. 
In one case, ice was removed more quickly and clearly 
from the treated than from the untreated surface after 
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Fic. 9. Ice removal from treated and untreated De-Icer sur- 
face on Mt. Washington. Icex VIa treatment. 
Icer. Treated surface on right. 


being exposed to over a week of normal flying weather. 
In another case the material was applied on a rainy day 
and flown through a heavy rainstorm for a number of 
hours without any apparent damage to the effective- 
ness as shown when ice was encountered. In other flight 
tests the treated portion of the De-Icers completely re- 
moved the accumulated ice at the end of one cycle, while 
additional inflations were required to remove the ice 
from the untreated area. 

It is felt that further trials may prove that one appli- 
cation will last through a winter season. In any event 
repeated applications of the silane or silicone can easily 
be made. 
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A Procedure for Determining the Critical 
Load for a Column of Varying Section 


G. H. MARTIN* 
Illinois Institute of Technology 


SUMMARY 


This manuscript presents an approximate method for deter- 
mining the critical load for a pinned end column of varying cross 
section. Although there are various procedures available for ob- 
taining a solution to the problem, the author believes that the 
application of the method suggested here is simpler in principle 
and that it requires relatively less calculation for the degree of 
accuracy obtained. 

The following examples of pinned end columns have been 
solved in the paper: 

(1) Columns of constant section. 

(2) Symmetrical columns of varying section. 

(3) Unsymmetrical columns of varying section. 


INTRODUCTION 


ie SLENDERNESS RATIO for a structural member 
that is subjected to axial compressive loads is de- 
fined as the ratio of L/k, where L is the length of the 
member and & is the least radius of gyration of its cross 
section. Members for which the value of L/k is greater 
than 160 may be called long columns. All columns re- 
ferred to in this paper will be considered as long col- 
umns. ‘ 

In Fig. 1 is shown a free end column subjected to an 
axial load P. The centerline of the column when in the 
deflected position is called the elastic curve, and the de- 
flection of any point on this curve from the straight line 
AB is designated as y. The shape of the elastic curve 
is governed by the familiar differential equation 


EI(d?y/dx®) = —M (1) 


In Eq. (1), E is the modulus of elasticity of the material, 
Tis the moment of inertia of the cross section, d*y/dx? 
is the rate of change in the slope of the elastic curve, 
and M is the bending moment at any point on the elas- 
tic curve. 

For the column shown in Fig. 1, the bending moment 
at a particular point on the elastic curve is Py. If this 
value for the bending moment is substituted for M in 
Eq. (1), we get 


EI (d*y/dx?) = —Py (2) 


The term critical load, as applied to an axially loaded 
column, is the load that will just cause buckling to be- 
gin. An extremely slight increase in the load above this 
value which causes a slight buckling, will increase the 
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deflection greatly. The fiber stress will increase with 
the deflection, and, when it has passed the yield point 
of the material, the column will fail with no further in- 
crease of the load. The critical load, therefore, is taken 
as the ultimate load and is represented by P in Eq. (2). 

If a column is considered as being composed of a 
number of sections and if the deflection curve is assumed 
to exist in the form of a straight line between sections, 
an expression may be obtained for the change in slope 
at any section. 

Let Fig. 2 represent a portion of a deflection curve, 
and consider the length of the column to be divided into 
a number of equal segments of length d. The slope of 
this approximate deflection curve midway between sec- 
tions a and b may be called y,’, then 
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Fic. 2. Portion of a deflection curve. 
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Ya’ = (M — Ya)/d 


Similarly the slope of the approximate deflection curve 
midway between sections } and c may be called y,’, then 


Yo’ = (Ye — My) /d 


These values for y,’ and y,’ may be considered as the 
average values of the slope of the actual deflection curve 
between sections a and b and sections b and c, respec- 
tively. The change in the average slope of the actual 
deflection curve as we move from section a’ to section 
b’ is 

prlauLey MOLEM i: Bete... BT 2% + Ya 

‘ ; ey d d 





Also, y” = (d/dx)(y’) and 
Ay’/ Ax = (y»" an: ya’) /d 


where d is the distance between sections a’ and b’. Sub- 
stituting the value of y,’ — y,’ into the last equation, 


Ay’/ Ax = (Ye — 2% + Ya)/d? 


This expression is an approximation to the average 
value of the second derivative which exists between 
sections a’ and b’. Assuming this as the value for 
*y/dx? at b on the actual deflection curve, Eq. (2) be- 
comes 


EI (2% — ye — Ya)/d?] = Py (3) 


The right-hand member of Eq. (3) represents the 
bending moment, and the left member represents the 
resisting moment at section 6. Under the critical load 
the resisting moment must equal the bending moment 
for equilibrium, and therefore Eq. (3) must hold for all 
sections along the column. 

Rewriting Eq. (3), 

(EI/ynd?)(2¥ — Ye — Va) = P (4) 


If we let P in Eq. (4) represent the critical load and 
since E and d are constants, then, for Eq. (4) to hold, 
the quantity 


(1/¥») (2% — Ye — Va) 
must be constant for all points along the column. Then 
(1/y») (2% — Ye — Ya) = K (5) 
The critical load; P,,, will then be 
Pa = K(E/d?’) 


COLUMN OF CONSTANT SECTION 


Before an attempt was made to apply the method to 
columns of varying section, it was tried on a column of 
constant section. The critical load for a column of ¢on- 
stant section may be readily found using Euler’s equa- 
tion, which gives the critical load as 


P. = 9.870EI/L? 
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Certainly, if the method is valid, it should agree with 
the Euler result. Thus the column of constant sec- 
tion was treated as a first attempt in verifying the 
method. 

The procedure for a column of constant section is 
the same as that which will be explained for a symmet- 
tical column of varying section. It will subsequently 
be shown that the critical load is 


Py = K,E/(L/S)? : (6) 
A comparison of the values for P,, as determined by 


Eq. (6) for several values of S,* with P., using Euler's 
equation, appears in Table 1. 





TABLE 1 
Comparison of P. with Exact Value 





Euler Equation 


Peo = 9.870EI/L? S=4 S = 6 S = 16 
Per = K,EI/(L/S)? 9.50EI/L? 10.08EI/L?2 9.80EI/L? 
Per cent error 3.75 2.13 0.80 





SYMMETRICAL COLUMN OF VARYING SECTION 


Knowing that the deflection curve for a symmetrical 
column is also symmetrical, we may consider just one- 
half of the column (see Fig. 3). 
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Fic. 3. Symmetrical column. 


The procedure is as follows: 

(1) Divide L/2 into a number of equal segments of 
length d. 

(2) Assume the deflection at the center to be 1.000, 
and assume the deflection curve to be a straight line 
from 0 to M as indicated by the dotted line in the figure. 

(3) Using Eq. (5) calculate K for each section except 
section 0, taking the deflections which correspond to 
the straight line OM. In applying the equation to a 
particular section, the value of the deflection at the 
section is the value to be used for y, in the equation. 
The value of the deflection at the adjacent section to 
the right is to be used for y,, and the value of the de- 
flection at the adjacent section to the left is to be used 
for y, in the equation. The value of K at section 0 is 
indeterminate, since the deflection at 0 is zero. 

(4) Average the K’s and let the value obtained be 
called Kg. 


* S is the number of equal segments considered over the full 
length of the column. 
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(5) Apply Eq. (5) to the center section M and then 
to each section to the left of M in succession, solving 
for y, in each case by using K,, as the value for K in 
Eq. (5). All along the column the value to be used for 
4 at a particular section is to be the value previously 
found for y, in the treatment of the immediate section 
to the right. 

(6) When the equation has been applied to the sec- 
tion adjacent to 0, the value determined for y, will be 
the deflection at point 0. Let the deflection at 0 be 
called yo. Since the load is applied at the ends where 
the deflection is zero, yo must come out equal to zero. 
A positive value for yo would indicate that K,, is too 
small. A negative value for yo would indicate that K,, 
is too large. The value of yo is affected in this way by 
K because a small value for K produces a small amount 
of curvature in the deflection curve; likewise, a large 
value for K produces a large amount of curvature in 
the deflection curve. The use of K,, as explained in 
step (5) will yield a value for yo near zero. For this 
reason K,, is a good value to use for K to start with. 
If several values of K are chosen and yp is calculated 
for each, K may be plotted against yo on coordinate 
paper. A smooth curve drawn through several plotted 
points will indicate the value of K which will make 
¥o equal zero (see Fig. 4). 





-Yo O t Yo 
Fic. 4. K plotted against yo. 








Let K, be the value of K corresponding to yo = 0. 
Then 


P., = K,E/(L/S)? 


where S is the number of segments in the full length of 
the column. The greater the number of segments used, 
the more accurate will be the value obtained for P,,. 

In Fig. 5 a symmetrical column with a constant 
moment of inertia I’ along its middle is shown. I is the 
value of the moment of inertia at the ends. /P,, has 
been determined for this type of column by the method 
just described, assuming the cross section to be circular, 
for a value of a/L equal to 0.4 and for a value of J/I’ 
equal to 0.6. The calculations appear in the Appendix. 

The critical load has been calculated for this column 
using various numbers of segments. Fig. 6 shows how 
the value of P,, approaches the exact value that has 
been determined by A. N. Dinnik! and which is de- 
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Fic. 5. Symmetrical column of varying section. 
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Fic. 6. Effect of increasing the number of segments in deter- 
mining the critical load. 


scribed by Timoshenko.? The curve shows K, plotted 
against S, the number of segments. The curve indi- 
cates that the greater the number of seyments consid- 
ered, the more accurate will be the value for the critical 
load. 


AN ALTERNATE METHOD FOR DETERMINING 
THE CONSTANT K;, 


In the procedure just presented, it was explained that, 
if too large a value for K is tried, the result will be a 
deflection curve that crosses the x axis. It can readily 
be seen that the constant K, in Eq. (6) is analogous to 
the constant m* in the Euler equation, which gives 


Pp = n*x*EI/L? 
where 1 is 1, 2, or 3, etc., depending on the order of the 
buckling load. 
The deflection curves for the first three buckling loads 


for a column of constant cross section with pinned ends 
are shown in Fig. 7. In Fig. 7 ” has the value 1, 2, and 
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Fic. 7. Buckling loads for a pinned end column. 
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3 for the deflection curves (a), (b), and (c), respectively. 
It is easily seen that the value of m,? may be obtained 
from m,* by multiplying m,? by the square of the ratio 
of the length between zero points in (6) and the full 
length of the column. That is, 


ng? = my2{(L/2)*/L4] = 4[(L*/4)/L4] = 1 
Similarly, 
n.? = n{(L/3)*/L4] = 9[(L*/9)/L‘] = 1 


The value of n? which gives the lowest buckling load 
is therefore obtained from any larger value of n? by 
multiplying the larger value by the square of the ratio 
of the distance between zero points to the full length of 
the column. This same method may be used as an al- 
ternate means of obtaining an approximate value for 
K, in Eq. (6). If in step (6) of the procedure for a 
symmetrical column it is found that yo comes out nega- 
tive, the value for K which yields this result can be 
greatly improved by multiplying its value by the square 
of the ratio of the distance between zero points and the 
full length of the column. This method for adjusting 
K may also be used when treating unsymmetrical 
columns. Such columns will now be discussed. 


UNSYMMETRICAL COLUMN OF VARYING SECTION 


In Fig. 8 a column of varying section which is un- 
symmetrical with respect to its centerline is shown. 
The procedure is similar to that for a column that is 
symmetrical, except that the deflection curve for the 
entire length of column must be considered. The pro- 
cedure is as follows: 
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Fic. 9. Deflection curve for an unsymmetrical column of vary- 
ing section. 
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Fic. 11. Column with abrupt changes in section. 


(1) Divide the length into a number of equal seg- 
ments of length d (see Fig. 9). 

(2) Let 0 designate the left end of the column and V 
the right end. If the deflection at N is taken as zero 
and the deflection at the section adjacent to N as 1.000, 
a deflection curve may be found which passes through 
these two points and also through point 0. 

(3) The critical load is given by Eq. (6) which gives 

P., = K,E/(L/S)? 

The approximate magnitude of K, may be found by as- 
suming a straight-line deflection curve from each end 
to the center temporarily (see Fig. 10). Assume the 
deflection at the center to be 1.000; then the deflections 
at the other sections can readily be determined from 
the straight-line relationship. Apply Eq. (5) to each 
section except sections 0 and N to determine K for each 
section except these two. 

(4) Average the values of K found in step (3) to 
obtain Kp. 

(5) Apply Eq. (5) to the section adjacent to N using 
Ka as the value for K in the equation. Solve for y,, 
taking y, as 1.000 and y, as zero (see Fig. 9). 

(6) Eq. (5) may then be applied to each of the re- 
maining sections to the left in succession, solving for 
Yq each time. The value used for y, at a particular sec- 
tion is to be the value found for y, in the treatment of 
the immediate section to the right. 

(7) When the section adjacent to 0 is reached, the 
value determined for y, will be the deflection at point 0. 
Let the deflection at 0 be called yo. Since yo must equal 
zero for the critical condition, a positive value for yo 
would indicate that K,, is too small and a negative 
value for yp would indicate that K,, is too large. The 
reason for this has been explained in step (6) in the 
procedure for a symmetrical column of varying section. 
Kg, will give a value of yo near zero. If several values 
of K are chosen and yy is calculated for each, K may be 
plotted against yo on coordinate paper just as it is 
shown in Fig. 4 for a symmetrical column. A smooth 
curve drawn through several plotted points will indi- 
cate the value of K which will make yp equal zero. Let 

K, be the value of K corresponding to yo equal zero. 
Then, 


PF sad K,E/(L/S)?* 


COLUMN WITH ABRUPT CHANGES IN SECTION 


The column with abrupt changes in section may be 
treated in the same manner as the column with gradual 
change in section. A column with abrupt changes in 
section is shown in Fig. 11. In dividing the length of the 
column into segments of equal length, it is advisable to 
take segments of such a length that sections within 
proximity of each discontinuity of cross section shall be 
considered. The larger the number of sections consid- 
ered, the greater will be the accuracy in determining the 
critical load. 
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Fic. 12. Symmetrical column of varying section. 


CONCLUSIONS 


The procedure presented in this manuscript for col- 
umns of varying section has been shown to be simple 
and to yield accurate results. The maximum error in 
determining the critical load for a variety of cases of 
symmetrical columns considered was found to be only 
1°*/190 per cent. This was the error when just four sec- 
tions were considered over the half length of the column. 

It has also been demonstrated that greater accuracy 
may be obtained in determining the critical load by con- 
sidering more sections. The effect of the number of 
sections considered on the per cent error in the critical 
load for a column of constant section is shown in Table 
9 


a 





TABLE 2 
Effect of Considering More Sections on the Accuracy of P-, 





No. of sections in L/2 2 3 5 8 
Per cent error in P., 3.75 2.13 1.12 0.80 





While the method has been applied to both symmetri- 
cal and unsymmetrical columns, in comparison to other 
existing methods it has been found to be especially ex- 
pedient for solving unsymmetrical columns, in that less 
calculation is required for an equivalent degree of ac- 
curacy in the critical load. 
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Fic. 13. K plotted against yo. 
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APPENDIX—COMPUTATIONS FOR DETERMINING THE 
CRITICAL LOAD FOR A SYMMETRICAL COLUMN OF.VARY- ° 
ING SECTION 


Circular section, a/L = 0.4 and J/I’ = 0.6. Let D 


be the diameter of the cross section. 


D, = 0.880 I, = 0.0294 
Dz = 0.930 I; = 0.0363 
D; = 0.980 I; = 0.0451 
D, = 1.000 I, = 0.0491 
D; = 1.000 I; = 0.0491 


Assuming a straight-line deflection curve from sec- 
tion 1 to section 5, y; = 0, ye = 0.250, vs = 0.500, v4 = 
0.750, and ys = 1.000. 


Ky = (0.0491/1.000)(2.000 — 0.750 — 0.750) = 0.0246 
K, = Ks = Ke = 0 
Ka = 0.0246/4 = 0.00615 


First calculation of deflections: 


K = (0.0491/1.000)(2.000 — y4 — 94) = 0.00615 


ys = 0.938 
K = (0.0491/0.938)(1.876 — 1.000 — ys) = 0.00615 
K = (0.0451/0:761)(1.522 — 0.938 — ye) = 0.00615 
MAL 0.480 

K = (0.0363/0.480) (0.960 — 0.761 — y1) = 0.00615 
yi = 0.118 
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Second calculation of deflections, trying K = 0.00690: 


K = (0.0491/1.000)(2.000 — +4 — ys) = 0.00690 
ys = 0.930 
K = (0.0491/0.930)(1.860 — 1.000 — ys) = 0.00690 
leas 0.731 
K = (0.0451/0.731)(1.462 — 0.930 — ye) = 0.00690 
ye = 0.420 
K = (0.0363/0.420)(0.840 — 0.731 — y:) = 0.00690 
y1 = 0.0280 


Third calculation of deflections, trying K = 0.00725: 
0 


725 


K = (0.0491/1.000)(2.000 — »,. — y.) = 0.0 


% = 0.925 
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K = (0.0491/0.925)(1.850 — 1.000 — ys) = 0.00725 

ys = 0.713 

K = (0.0451/0.713)(1.426 — 0.925 — ye) = 0.00725 

ye = 0.386 

K = (0.0363/0.386)(0.772 — 0.713 — y1) = 0.00725 

v¥, = —0.0181 

= K,E/(L/S)?° 

0.00718 E/(L/8)?, for I’ = 0.0491 

0.00718 EI'/0.0491(L/8)? = 9.35EI'/L?, for any 
value of I’ 


Ny 
I | 


From Dinnik’s equation, P,, = 9.45EI’/L?. 
Per cent error: 


(9.45 — 9.35)/9.45]100 = 1.06 


Letters to the Editor 


Dear Sir: 

A “‘beat”’ occurs when two simple harmonic motions, or two 
sine functions of any amplitudes but only slightly different fre- 
quencies, are added together. The result is a wave that is ap- 
proximately sinusoidal, except that its amplitude and frequency 
vary periodically. It is well known how to find the beat fre- 
quency, wz, of the amplitude (and frequency) variation and the 
amplitudes of the two components from the composite curve. 
This letter shows that the frequencies of the components may also 
be found from the curve as follows: The frequency of the com- 
ponent with the larger amplitude, w;, is the same as the average 
frequency of the composite when measured over an integral 
number of beat cycles. The frequency we of the smaller com- 
ponent is w, + we if the apparent frequency of the composite 
curve is larger at the portion of the beat cycle having a large 
amplitude than it is at the portion having a small amplitude 
(i.e., if the peaks are spaced closer together at the wide part of the 
beat than at the narrow part). w is w, — wa if the peaks are 
spaced farther apart at the wide part of the beat than at the 
narrow part. 

Proof: To prove the stated properties, we study the synthesis of 
a beat, adding the two component sine waves: 


x = A, COS at 
Ag COs (wet + 4) 


Xp = xX; + Xe 


X2 


where A, is the larger amplitude and @ is the phase lead of the 
smaller component relative to the larger at time ¢ = 0. Fig. 1 
gives the addition of the two sine waves in a vector diagram at 
an arbitrary instant. Each vector rotates counterclockwise with 
its own frequency, and its projection on the X axis is x. The 
larger vector A, is drawn from the origin. The two component 
vectors go around approximately together, changing their relative 
phase angle @ by only a few degrees in each revolution. Hence, 
the amplitude of the sum stays fairly constant, and in each cycle 
the composite motion closely resembles a sine wave. 
The phase angle between the two components is 


8 = (wot + %)~% = (we — wy)t + % = wat + % 


we is the beat frequency and the frequency with which the two 
vectors come in phase. 

Fig. 2 shows the variation of the resultant vector B relative to 
the larger component as @ changes with time. Fig. 3 gives the 
composite curve. 


The “‘envelope”’ of the beat will be defined as the length of the 
vector B (Fig. 2) versust. (This also happens to be the smooth- 
est curve tangent to the composite curve.) We can find the 
amplitudes of the components by measuring 2B,,,,. and 2B,,;,. 
on the envelope. Then, 


A, eat 7 | ae + 2Brmin,)/4 
and 
A, = (2B mez. ba 2Brmin.)/4 


The frequency of the larger component may be found by meas- 
uring the frequency of the composite wave averaged over one 
cycle, from maximum to maximum of the beat, because the com- 
posite amplitude vector B never gets more than 90° out of phase 
with the larger component vector A;, and B returns to the same 
phase relative to A, after each beat cycle (see Fig. 2). The fre- 
quency w; is then the number of radians of rotation (w)t2 — wyt;) 
of the vectors divided by the corresponding length of time (see 
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Fic. 1 (left). Vector addition of components A; and A» giving 
beat vector B. Fic. 2. (right). Variation of resultant vector B 
relative to the larger component A). 
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Fic. 3. Record of beat with frequency calculations. 
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LETTERS TO 


example, Fig. 3). It will usually be sufficiently accurate to 
simplify .the calculation by taking the largest peaks as the limits 
of the interval and calling wt; = 0 and wf. equal to an even mul- 
tiple of 360° or 27 radians. 

Now, to find the frequency of the smaller amplitude component, 
we note that, when it has a higher frequency than the larger 
component, A: is rotating faster than A, (in Fig. 1), and the com- 
posite vector B must be rotating faster than A, when B is large 
and slower than A, when B is small. In other words, when the 
peaks are more closely spaced and have a shorter period at the 
maximum width than they do at the minimum width, as in Fig. 3, 
we know that w. = w, + wg. Similarly, if the peaks are less 
closely spaced at the maximum width than at the minimum, 
the composite vector B must be rotating more slowly than A, 
when it is large and faster than A; when B is small, and Az must 
be rotating more slowly than A;. Then w. = w, — we if we is 
considered to be a positive number. 

Applications: The writers developed the analysis when numer- 
ous problems involving beats presented themselves during the 
analysis of aircraft vibrations. The more important cases were 
separating response to a vibrator from ‘‘hash”’ response to engine 
or propeller excitation near the vibrator frequency and separating 
normal modes in records of free vibration. In the latter case 
the measured composite curve is modified by the multiplication 
by an exponential function of time as the motion decays due to 
damping; the record could be made into a beat between constant 
amplitude components by introducing automatic volume control 
if the damping were the same in the various modes. 

C. D. P. Pancu 

Curtiss-Wright Research Laboratory 
CHARLES C. KENNEDY 

Howard University 





Dear Sir: 

In the paper by J. G. Sawyer, ‘‘Design and Application of 
High-Pressure Coefficient Axial-Flow Fans’’ (JOURNAL OF THE 
AERONAUTICAL SCIENCES, Vol. 12, No. 3, p. 281, July, 1945), 
the following points seem to deserve attention: 

(1) Mr. Sawyer’s introductory remarks that we now have 
available design methods that have made possible axial-flow 
fans with pressure characteristics superior to all other types of 
fans is certainly not true, unless Mr. Sawyer is thinking of the 
other fans as being of the centrifugal type, in which case the 
writer agrees with him. However, exhaustive investigation 
undertaken during recent years by the writer in collaboration 
with V. H. Pavlecka, has demonstrated that a certain new type 
of compressor has superior pressure and also pressure-volume 
characteristics to the axial-flow, as well as to the centrifugal- 
flow, compressor. ; 

(2) The writer objects to calling AP the fan pressure rise, 
since, in reality, because of unavoidable friction losses, the true 
pressure rise is somewhat less than that given by Eq. (3). The 
true pressure rise can always be obtained by multiplying Eq. (3) 
by the internal efficiency of the fan. 

(3) The writer furthermore objects to the use of inconsistent 
nomenclature employed in the Eqs. (1) to (4). 

(4) Sawyer’s Eq. (2) is not a general equation for the pressure 
tise; the general equation for the pressure rise in terms of-the 
airfoil lift coefficient is 


(CLWeolpw2/4r)(Cm/u) cot Bo 


where sin Bo = Cn/Wa. It is apparent that Sawyer’s Eq. (2) 
represents the special case in which 


(Cm/u) cot Bo = 1.0 
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This expression establishes a fixed relation between the velocity 
components of the velocity vector diagram, viz.: 


Cu aoa — Cus 


Sawyer does not specify that Eq. (2) holds only for the special 
case in which Cy, = —C,,. Sawyer’s Eq. (2) is apparently 
taken from the book by Curt Keller, Eq. (28), page 9, the latter 
being correct for the conditions that Keller carefully points out, 
viz., Cu = 0 and AP, = stage pressure rise—i.e., the total 
pressure rise over the guide vanes and the fan wheel as well. 

(5) Sawyer’s Eq. (4) is not the general equation for the fan 
pressure rise but is true only for the special case in which Cy, = 
—Cy2, as in Eq. (2). The general equation for the fan pressure 
rise is 

AP = puACy — (p/2) AC, — pCuACy 


Sawyer’s statement, therefore, that the only means of increasing 
AP is to increase AC, is certainly not true, since AP may be 
increased or decreased by decreasing or increasing C, as well as 
AC,. In fact, AP may even decrease when AC, alone is 
increased. This will occur when AC, > (u — Cy). It is, there- 
fore, not necessarily true that higher pressures call for highly 
curved blades, since higher pressures may be obtained by de- 
creasing C,y,:—i.e., by allowing the blades to operate at smaller 
pitch angles. Thus, flat blades may develop high pressures, as 
has been done in a number of machines built and tested in 
Europe. 

(6) According to Fig. 2, Cy. = 0 and the correct equation 
for the fan pressure rise becomes 


AP = puAC, — (p/2) AC,? 


Then, using Sawyer’s definition for the pressure coefficient, the 
correct expression for the pressure coefficient, when Cy, = 0, is 


vy = 2(AC,/u) — (AC,/u)? 


The maximum pressure coefficient obtainable is y = 1.0. This 
result may be obtained by the method generally employed in 
the elementary calculus for finding the maximum value of a 
function. The fans that Sawyer claims have developed pressure 
coefficients greater than 1.0 certainly did not operate with 
Cu. = 0 as indicated in Fig. 2. 

(7) For the special case when 


AP = pudAC, 


it is possible to obtain pressure coefficients greater than 1.0, 
although the writer doubts that values of y = 3.0, as indicated 
on Fig. 7, can be obtained with the airfoil sections described 
and shown by Sawyer. The following example may help to 
clarify this statement: 

Example: For maximum efficiency, the aerodynamic theory 
of axial-fan design teaches that 8. should be equal to 45°. In 
the interest of efficiency it is assumed that Bo = 45°. The 
following expressions for the relative inlet and outlet angles can 
be derived: 


tan 8; = 1/[1 + (¥/4)] 


and 

tan 6, = 1/[1 — (¥/4)] 
For yy = 3.0, 8; = 29° 45’ and 8, = 75°58’. The turning angle 
is 

v = 6 — B, = 46°13’ 


It has been the writér’s experience that turning angles 3 greater 
than 20° could not be attained because of flow separation so that 
values of 8 = 46° 13’ appear to be outside the range of actual 
attainment by known means in simple cascades. 
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(8) In Fig. 3 of the author’s paper, the velocity vector dia- 
grams are incorrectly drawn. Fig. 3 shows that Cui = 0 and 
Cuz = 0. Obviously, therefore, the fan cannot absorb any 
power (since AC, = 0) and, hence, will produce no flow, since 
it will not rotate. In the referenced N.A.C.A. publication 
the velocity vector diagrams for the Schicht fan are shown 
correctly. 

(9) In Fig. 7 Sawyer attempts to show the variation of 
pressure with flow for various values of y. What is actually 
shown is certainly not clear since all curves for y = constant 
will have a positive slope and not any slope as is indicated. This 
fact can be shown by the following analysis: 

The pressure rise may be expressed by the relationship, Eq. (1) 


AP = (p/2)u*y 


Upon differentiating and dividing both sides of the resulting 
equation by the differential change of volume dV, we get 





ap duo, Wy 
cae” 6s ae 


New for ¥ = constant, the above equation reduces to 
dAP/dV = puy(du/dV) 


This equation shows that for positive values of y the slope of the 
pressure-volume curve dAP/dV is always greater than zero, 
since for a constant value of y the characteristics of ail fans are 
such that du/dV > 0—i.e., when the speed increases the volume 
also increases. All curves shown in Fig. 7 should therefore be 
drawn with a positive slope. 

It is the general practice to represent fan characteristics at 
constant speed. Thus, for constant speed we have 


dAP/dV = (p/2)u*(dy/dV) 
The volumetric flow V may be expressed by the relationship 
V = Fou 


where F = flow area, 


@ = C,,/u = flow coefficient 
Hence, for u = constant 
dAP/dV = (pu/2F)(dy/d¢) 


‘Most fans show a large negative slope of dy/d@ at the design 
value of ¢. However, there is a functional relationship between 
y and ¢ which gives dy/d@ > 0 at the design value of ¢. When 
Cu. < 0, it can be shown that dy/d¢ > 0. However, it is im- 
portant to note that even with Cy, < 0, dy/dd may be negative 
provided C,, is not sufficiently great. This is due to the friction 
losses that Sawyer chooses to neglect but which, unfortunately, 
are always present in practice. 

Sawyer fails to make clear this well-known fan characteristic, 
and the presentation of Fig. 7 with the incorrectly shown nega- 
tive slopes seems only to add to the confusion. 

(10) Sawyer correctly indicates and appreciates the design 
principles as pointed out for the first time in print by Sérensen— 
viz., that high pressure fans or compressors can be designed 
rationally only on the basis of the channel theory familiar from 
the turbine design, instead of on the basis of the airfoil theory of 
lift. The writer was active in the development of axial com- 
pressors before the appearance of Sérensen’s paper in this country 
and recalls the difficulties he had in convincing several eminent 
aerodynamicists to accept the channel theory in lieu of the airfoil 
theory for compressor design, a conclusion he had arrived at 
independently by his own analytical work at that time. Because 
of the lack of understanding of the subject and insufficient 
imagination on the part of aerodynamicists into whose hands the 
development of the axial compressor has been unfortunately 
entrusted and the confusion contributed by the modern aero- 
dynamic theory by veiling the true nature of the art of compres- 
sion of gases, the developmental progress of the axial compres- 
sors has been substantially delayed in this country—a dis- 
advantage the steam turbine luckily did not suffer from during 
its formative development some 40 years.ago. 

FREDERICK DALLENBACH 
8918 7th Ave. 
Inglewood, Calif. 
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Aero Medical Aspects of Cabin Pressuriza- 
tion for Military and Commercial Aircraft 


W. R. LOVELACE, II,* ann A. P. GAGGEt 
Wright Field 


SUMMARY 

Whenever aircraft fly at high altitudes for long periods of time 
and over long distances, the foremost design requirement for air- 
crew efficiency and passenger comfort is cabin pressurization. In 
military aircraft, cabin pressurization is used to eliminate com- 
pletely the need for continuous use of oxygen, to prevent aero- 
embolism only, or simply to reduce the cabin to a level where use 
of pressure breathing is not required. In commercial aircraft, the 
health, comfort, and safety of the passengers are of prime con- 
sideration. In military aircraft, a significant hazard introduced 
by cabin pressurization is explosive decompression caused by 
structural failure of the cabin when under enemy gunfire. Criteria 
are presented by which safe operating differentials may be chosen 
to minimize this hazard. The advantages of isobaric cabin al- 
titudes below the now conventional 8,000-ft. level for military 
aircraft are reviewed as they apply to passengers of commercial 
aircraft. Finally, emergency procedures necessary for crew and 
passenger safety in case of failure of the cabin-pressurization sys- 
tem are discussed. 


INTRODUCTION 


2 apras THE PRESENT WAR, technical advances toward 
realization of high-altitude flight have been many. 
The three most important on the engineering side 
are the development of the turbosupercharger and jet 
propulsion engines and the evolution of the laminar 
flow wing. Now within reach are flights by military 
aircraft well above limits tolerated by human beings 
breathing pure oxygen. Sustained flights at high al- 
titude are also possible for both military and commer- 
cial aircraft. Parallel with these rapid engineering 
advances, methods for maintaining air-crew efficiency 
and safety as well as passenger comfort must be con- 
tinually developed. The most logical step taken by 
the aircraft designers in this direction has been cabin 
pressurization. It is the purpose of the present paper 
to present the aero medical requirements for health, 
comfort, efficiency, and safety of air crews and passen- 
gers while flying in pressurized-cabin aircraft. 


PHYSIOLOGIC REQUIREMENTS IN FLIGHT 


It is a well-accepted fact, based on long flying and 
operational experience, that optimal efficiency and 


Presented at the Los Angeles Meeting, I.A.S., Los Angeles, 
August 16-17, 1945. 

* Colonel, Medical Corps. Chief, Aero Medical Laboratory, 
Air Technical Service Command. On leave of absence from the 
Mayo Clinic. 

t Lieutenant Colonel, Army Corps. Chief, Biophysics Branch, 
Aero Medical Laboratory, Air Technical Service Command. 
Formerly at the John B. Pierce Laboratory, Yale University. 


safety of air crews require the supplementary use of 
oxygen above 10,000 ft. for daytime flight and above 
5,000 ft. for night flight. Flight without oxygen above 
these levels will result in increasing degrees of anoxia 
and in progressive loss of judgment and flying efficiency. 
The demand oxygen system now used by the Army Air 
Forces is designed to increase automatically the con- 
centration of oxygen from that of pure air at 5,000 ft. 
to 100 per cent oxygen at 33,000 ft., so that the concen- 
tration in the lungs remains equivalent to that at 5,000 
ft. or below. Above 33,000 ft. the margin of safety 
using pure oxygen decreases steadily to 40,000 ft., 
which is the oxygen altitude equivalent to 10,000 ft. 
with air. Above 40,000 ft. the 10,000-ft. equivalent 
can be maintained only by the use of pressure breathing. 
The top level now recommended for continuous use of 
pressure breathing is 42,000 ft. with 6'/2 in. of water 
positive pressure; this combination corresponds to the 
10,000-ft. air equivalent. Under emergency conditions 
for short periods, pressure breathing may be used to 
50,000 ft. with 12 in. of water positive pressure, corre- 
sponding to 16,000-ft. air equivalent. For emergencies 
in flight above 50,000 ft., a lightweight quickly inflatable - 
pressure suit would be required to reduce the equiva- 
lent altitude about the flier to 40,000 ft. or below. 

Statements similar to those presented in the preceding 
paragraph apply to the continuous flow type oxygen 
system. However, where the demand system supplies 
the necessary supplementary oxygen for all levels of 
exercise, the level of flow in the former system must be 
increased as the level of exercise increases in order to 
maintain a fixed safety level. At altitudes above 
25,000 ft. and for the same degree of safety, the de- 
mand system is much more economical in the use of 
oxygen than the continuous flow. 

Aero-embolism, which manifests itself in the form of 
bends, of chokes, and of skin and neurologic symptoms, 
rarely occurs below 25,000 ft. However, it becomes in- 
creasingly significant as the altitude increases above 
30,000 ft. For 1 hour’s exposure or less at 35,000 ft., 
one person in ten would be likely to be incapacitated; 
one in four, at 40,000 ft. Few persons can stay more 
than 20 min. above 40,000 ft. without suffering from 
some form of aero-embolism. The incidence of aero- 
embolism is increased by exercise on the part of air 
crews. The incidence of aero-embolism can be greatly 
reduced by prebreathing 100 per cent oxygen for '/: 
hour before flight and breathing 100 per cent oxygen 
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for the entire flight. There is a small percentage of 
flying personnel for whom oxygen prebreathing gives 
no protection from aero-embolism, who suffer regularly 
on flight above 35,000 ft., and who should be disquali- 
fied for high-altitude flight. 

At altitudes below 10,000 ft., proficiency in night 
vision and dark adaptation is an extremely important 
consideration, since it is greatly affected by mild de- 
grees of anoxia. For example, to maintain an equiva- 
lent stimulus on, and response from, the retina, the 
percentage of increase in light intensity for various al- 
titudes is as follows: 


Sea level O per cent 

5,000 ft. 23 per cent 
10,000 ft. 59 per cent 
12,000 ft. 78 per cent 
14,000 ft. 101 per cent 
16,000 ft. 140 per cent 


Although use of high concentrations of oxygen from 
sea level up is desirable for the preservation of night 
vision, the small loss of night visual efficiency at 5,000 
ft. is the maximum allowable for effective operational 
purposes. Above 5,000 ft., use of oxygen is a definite 
requirement to preserve night vision. A second visual 
phenomenon associated with mild anoxia is the dim- 
ming of peripheral vision, especially when illumination 
is poor. Both Hecht and Livingston have shown 
that this effect becomes increasingly serious above 
6,000 ft. 


LEVELS OF CABIN PRESSURIZATION 


There are essentially three types of human design 
requirements for use of cabin pressurization: 

First, the cabin shall be maintained at levels where 
continuous use of oxygen equipment will not be needed. 

Second, when oxygen is used continuously, the cabin 
shall be maintained at a level low enough to prevent 
incidence of aero-embolism. 

Third, for flights above 40,000 ft. , sufficient pressuri- 
zation shall be used to obviate the use of pressure 
breathing and the risk of mild anoxia above 42,000 ft. 
with pressure breathing. 

To meet these human requirements the engineer 
must stress his cabin fuselage to take the necessary 
differential pressure. The maximal required differential 
pressure is set by the operational ceiling of the aircraft 
(in contradistinction to the service ceiling). 

The first of the foregoing requirements is used by the 
Army Air Forces in the design of its pressurized bomb- 
ers. For example, the B-29 is designed for a 6.55 Ibs. 
per sq.in. maximal operating differential, which allows 
an 8,000-ft. cabin altitude at 30,000 ft. or a 10,000-ft. 
cabin altitude at 35,000 ft. 

The second of the foregoing requirements is followed 
in the design of pressurized fighters. The differential 
pressure most frequently chosen by the aircraft design- 
ers is 2.75 Ibs. per sq.in., which allows a 25,000-ft. 
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cabin altitude at 40,000 ft. or a 30,000-ft. cabin altitude 
at 50,000 ft. 

Choice of any design cabin differentials between the 
2.75 Ibs. per sq.in. and the 6.55 Ibs. per sq.in. value, 
the aircraft engineer will quickly note, does not increase 
the tactical efficiency of air crews when flights above 
30,000 ft. are possible. A value of 4.0 lbs. per sq.in. at 
30,000 ft. would result in a 15,000-ft. cabin altitude, at 
which level oxygen would have to be used continuously. 
At 40,000 ft. a 2.75 lbs. per sq.in. differential is adequate 
for the second requirement; the possible saving in 
oxygen consumption using a 4.0 Ibs. per sq.in. differen- 
tial instead of a 2.75 Ibs. per sq.in. would be balanced 
by the increased loading of the cabin superchargers and 
heavier construction of the cabin walls and pressure 
seals. 

The third requirement, listed previously, would be 
used only under emergency or combat conditions, when 
other hazards, such as explosive decompression (which 
will be discussed later), prevent the use of either of the 
first two requirements. 

In Fig. 1, the various factors presented in the two 
previous sections are outlined to indicate their interre- 
lationship. 


EXPLOSIVE DECOMPRESSION 


The principal hazard to air crews of tactical aircraft, 
when operating with pressurized cabins, is the possible 
damage caused by the expansion of internal body gases 
as a result of explosive or rapid decompression of the 
cabin itself. The failure of the cabin structure may 
have been caused by enemy gunfire. The degree of 
explosive decompression that one can withstand safely 
is determined by the extent and the rate of expansion of 
the internal body gases such as are found in the stom- 
ach, intestines, and lungs. The experimental studies of 
Major H. M. Sweeney, of the Aero Medical Laboratory, 
Wright Field, summarizéd in detail in the October, 
1944, issue of the Air Surgeon's Bulletin, have shown 
that the two most significant factors associated with 
sudden decompression are the time of decompression 
(t) and the relative expansion of internal gases (RGE). 
The latter is evaluated by the following relation: 


RGE = (P, — 0.91)/(P. — 0.91) (1) 
where 
P. = cabin pressure before decompression in lbs. 
per sq.in. 
P, = ambient or final pressure after decompression 


in lbs. per sq.in. 
0.91 = vapor pressure of water at body temperature, 
98.6°F. in Ibs. per sq.in. 


In particular, it was demonstrated in the foregoing 
series of experiments that for rapid or instantaneous 
decompressions (approximately 0.01 sec. or less) a rela- 
tive gas expansion of 2.3 approached the maximum 
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tolerable. As the time of decompression increases, the 
increased pressure of expanding gases has time to be 
absorbed partially by neighboring tissues and the gases 
can escape from the body; thus the maximal tolerable 
RGE would increase. In Fig. 2, the maximal tolerable 
decompression conditions, observed for human subjects, 
have been plotted in terms of RGE and the time of de- 
compression, ¢. A straight line, drawn between the 
topmost series of experiments, is described by the 
equation: 


RGEmaz, = 2.1 + 17.08 (2) 


This line defines the borderline decompression condi- 
tions between the zone of safety, determined by ade- 
quate laboratory experience, and a zone of uncertainty, 
in which experimental data either are limited or are 
unknown. The line described by Eq. (2) can be consid- 
ered by the aircraft engineer a maximal design limit for 
tolerable decompression that includes a margin of safety 
sufficient for military operations. It should be added 
at this point that it is probable that certain conditions 
above this line are tolerable; however, the present 
limit is proposed until new data are available for which a 
new design line can be drawn. 

From Fliegner’s equations for the flow of air from an 
orifice, it can be shown that the time of decompression, 
t, in seconds, may be described with sufficient accuracy 
by the equation: 


t = 0.22(V./A) V(P, — P.)/Pa (3) 
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Fic. 2. Relation between the maximal allowable relative expan- 
sion of internal gases (RGE) and time of decompression. 


where 
V. = volume of pressurized cabin, cu.ft. 
A = cross-sectional area of hole caused by struc- 


tural failure, sq.in. 
P, and P, as in Eq. (1) 


A preliminary test of this equation by an actual de- 
compression in flight indicates that, because of a pos- 
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sible orifice coefficient less than unity, the observed 
times of decompression may be longer than calculated 
in Eq. (3). For the present analysis an equation giving 
the probable minimal decompression time for a given 
set of design factors is desirable. 

Combining Eqs. (2) and (3) the maximal tolerable 
relative expansion of internal gases is given by the re- 
lation: 


RGEmaz. = 2.1+3.8(V./A) V (P. — P)/P. (4) 


The criterion of safety for an air crew during any de- 
compression in a pressurized-cabin aircraft is that the 
RGEnaz, aS calculated by Eq. (4) is greater than or 
equal to RGE as calculated by Eq. (1). When RGE 
calculated by Eq. (1) exceeds RGEnar,, dangerous or 
questionable conditions exist. These criteria for safety 
should be applied at the maximal service ceiling of the 
aircraft under consideration. Under these circum- 
stances it can be demonstrated that flight conditions 
below the service ceiling will also satisfy the criteria 
with a further margin of safety. 

The factor A, the area of the fuselage surface lost by 
structural failure caused by enemy action, must in the 
long run be evaluated by gunfire tests. For preliminary 
design purposes it can be assumed to be the area of a 
jettisonable bubble canopy, a nose section, or an escape 
hatch door. Practical application of the foregoing 
criteria shows that: 

(a) For pressurized fighter aircraft where the cross- 
sectional area of the bubble canopy (approximately 
600 sq.in.) is large in relation to the cockpit volume 
(50 cu.ft.) the limit of the zone of safety is described by 
a choice of cabin altitudes or differential cabin pressures 
so that for any flight altitude the RGE does not exceed 
2.3. 

(b) For heavy bombardment aircraft (2,000 cu.ft. 
and up) flying with a cabin differential of 6.55 Ibs. per 
sq.in. at 40,000 ft., the largest allowable explosive ori- 
fice A is 7,000 sq.in. or a hole 8 ft. in diameter. Since 
this hole is ten times larger than that normally expected 
from enemy action, decompression may, as a rule, be 
considered a minor hazard for this type aircraft. 

(c) The greatest variability in allowable operating 
conditions occurs in medium bombardment aircraft 
with cabin volumes in the range of 200 to 400 cut. 
This is illustrated in the following typical example for 
an aircraft with a cabin volume of 350 cu.ft. flying at 
45,000 ft. 

(1) For loss of a jettisonable bubble canopy (2,000 
sq.in.), RGEmaz, is 2.75 and maximal allowable cabin 
differential is 2.15 lbs. per sq.in. 

(2) For loss of a nose section (1,000 sq.in.) RGE naz. 
is 3.85 and allowable differential is 3.5 Ibs. per sq.in. 

(3) For loss of an escape hatch (600 sq.in.) RGEnaz. 
is 5.7 and the allowable differential is 5.7 ibs. per 
sq.in 
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In Fig. 1, the line for an RGE of 2.3 applicable to 
fighter aircraft has been drawn to indicate the limit of 
safe conditions for explosive decompression. Exami- 
nation of this figure will show that for fighter aircraft 
in combat above 50,000 ft. the choice has to be made be- 
tween the dangers of explosive decompression and the 
increasing degrees of anoxia while using pressure breath- 
ing oxygen. 


AURAL AND SINUS PHENOMENA 


From decompression experiments on human subjects 
it has been observed that the pressure in the middle ear 
is relieved automatically for rapid changes in outside 
pressure. A normal ear can easily adjust itself 
to a rate of pressure loss of 100 Ibs. per sq.in. per see, 
without any aural or pressure sensation at all to the 
subject. This statement is not true for air-crew mem- 
bers with colds and respiratory infections. A passen- 
ger with an occluded sinus would suffer severely from 
pain when either the volume or the pressure of the 
trapped gases increased as little as 30 per cent, regard- 
less of the rate of decompression. In the design of 
cabin-pressure regulators, the Army Air Forces require 
that the drop in cabin pressure be limited to 1 Ib. per 
sq.in. per sec., a value that can be considered conserva- 
tive in the light of explosive decompression data now 
available. 

With a rise in cabin pressure, the negative pressure in 
the middle ear can be relieved only by a voluntary ef- 
fort on the part of the air crew. The ability to clear 
ears varies greatly with the individual and with his 
flying experience. For military aircraft the rate of rise 
of cabin pressure is limited by pressure control to 1 Ib. 
per sq.in. per min. (approximately 2,500 ft. per min.), 
provided the rate of dive of the aircraft itself does not 
exceed this value. For commercial aircraft, past expe- 
rience hes shown that,.0.10 lb. per sq.in. per min. or 
300 ft. per min. at sea level is the most comfortable de- 
scent rate for the average passenger. Recently, it was 
shown at the School of Aviation Medicine, Randolph 
Field, Tex., that, when the cabin is fluctuated between 
7,700 and 8,300 ft. at a.rate of 600 ft. per min. or 0.25 
Ib. per sq.in. per min., sensation on the ear drums was 
so small as not even to disturb sleep. This result will 
indicate the required constancy in the desired con- 
trols. 

When the cabin controls are fitted with a rate of climb 
and descent adjustment, discomfort to the passenger 
can be reduced to the minimum. This consideration is 
especially important to engineers designing aircraft for 
commercial air lines, who must design their pressure 
cabin controls so that passengers with upper respiratory 
infections and sinusitis will be comfortable in flight. 
With such controls it would be possible to maintain in 
the cabin a rate of descent lower than that of the plane 
itself, and thus an actual rapid descent to an airport is 
minimized. 





dec¢ 
to b 
enti 
cu.f 
expt 
that 
deb: 
Dur 
high 
the 


cons 
brok 
gunt 
stati 
W 
tial 
thos 
tact: 
be < 
are: 
(a 
shou 
port 
(b 
nets 
seng 
boar 
A 
rienc 
in pi 
shou 
layer 
lami: 
sing] 
erati 
Plex: 
of ex 
tion. 
Th 
the ¢ 
aircr: 
visua 
not 
Whe: 
mal t 
cepti 
expec 
optic 
in pa 
is re 
neces 
impr 
Ide 
cal ai 





sable to 
limit of 
Exami- 
aircraft 
1ade be- 
and the 
breath- 


ubjects 
idle ear 
outside 

itself 
er sec, 
to the 
’ mem- 
assen- 
y from 
of the 
egard- 
sign of 
‘equire 
lb. per 
iserva- 
a now 


sure in 
iry ef- 
» clear 
th his 
of rise 
> 1 Ib. 
min.), 
2s not 
expe- 
in. OF 
le de- 
t was 
dolph 
‘ween 
- 0.25 
S was 
t will 
con- 


‘limb 
nger 
on is 
t for 
ssure 
tory 
ight. 
in in 
lane 
rt is 





MEDICAL ASPECTS OF CABIN PRESSURIZATION 147 


SAFETY DESIGN OF PRESSURE-CABIN AIRCRAFT 


When a side sighting blister is lost in the B-29, the 
decompression time has been determined by firing tests 
to be approximately 1 sec. for the full operating differ- 
ential of 6.55 Ibs. per sq.in. and cabin volume of 2,200 
cu.ft. During this second, approximately 3,500 hp. is 
expended. From the safety point of view it is essential 
that this energy does not reappear in the form of flying 
debris, such as loose items of personal equipment. 
During a decompression in the B-29, air velocities as 
high as 140 m.p.h. have been observed momentarily in 
the interconnecting tunnel. A flier in the tunnel at 
the moment of decompression would be ejected with 
considerable force. Also, the rush of air through the 
broken blister has sufficient directed energy to blow a 
gunner overboard if he is not properly secured to his 
station by a safety belt. 

With proper indoctrination and training in the poten- 
tial hazards of cabin pressurization, accidents such as 
those described previously for the B-29 are of negligible 
tactical significance. However, the lessons that can 
be applied to passenger aircraft are obvious. These 
are: 

(a) Restricted passageways between compartments 
should be avoided. This precaution is especially im- 
portant as the size of the aircraft increases. 

(b) Regular exits should have safety harnesses or 
nets across the inner side during flight to prevent pas- 
sengers and crews from being accidentally blown over- 
board. 

Another important lesson from the operational expe- 
riences of the B-29 is that the use of single-layer Plexiglas 
in pressurized aircraft is an undesirable practice and 
should be avoided wherever possible. Since the single- 
layer Plexiglas blisters of the B-29 were replaced by the 
laminated type about a year ago, there has not been a 
single report of a blister failure of this type from op- 
erational causes. In addition to the use of laminated 
Plexiglas, another desirable practice is to avoid the use 
of extremely large Plexiglas sections without segmenta- 
tion. 

The increasing tendency to use formed Plexiglas in 
the design of nose sections and canopies of pressurized 
aircraft leads to a word of precaution. Until now, the 
visual properties of Plexiglas canopies and sections have 
not been so good as those of grade B plate glass. 
Where the line of sight is greater than 55° to the nor- 
mal to the Plexiglas surface, serious errors in depth per- 
ception and loss of acuity and of distant vision can be 
expected. Use of flat laminated plate glass of grade A 
optical quality is a requirement for crew compartments 
in passenger aircraft. For tactical aircraft where speed 
is required, a compromise, using curved plastic, is 
necessary. However, the plastic industries can well 
improve the visual properties of viewing sections. 

Ideally, all exit hatches should open outward in tacti- 
cal aircraft. At the present time this is true for fighter 


aircraft. When space for the air crews is restricted— 
as it would be, for example, in medium bombardment 
aircraft of 300 to 400-cu.ft. volume—outward opening 
of hatches is also desirable. For the heavy bomber air- 
craft, such as the B-29, and for passenger aircraft, 
where cabin space is not at such a premium, inward 
opening doors are allowable. A small residual cabin 
pressure can preclude the use of inward opening doors 
in an emergency. In actual practice, all inward open- 
ing emergency doors or exits should be able to be opened 
by simple mechanical means against at least a 0.5 Ib. 
per sq.in. inside pressure. 

Escape from high-speed aircraft is becoming of in- 
creasing concern to the aircraft designer and to flying 
personnel. Recently, it has been revealed that the 
German Air Force has developed an ejection seat pro- 
pelled by a small gunpowder charge expanding into a 
piston. This type of seat at the present time appears to 
be the most practical approach to an escape method 
from high-speed aircraft. When an ejection seat is 
used in pressurized-cabin aircraft, the canopy must be 
ejected simultaneously with the seat causing the pilot 
or air crew to experience a rapid or even explosive de- 
compression. It is therefore important that a differen- 
tial pressure in the cabin be so selected by an RGE con- 
trol that the safety criteria for decompression, as out- 
lined previously, always apply. In particular for 
fighter aircraft, for flight altitudes above 38,000 ft. 
(Fig. 1), an RGE control of 2.3 should override the 2.75 
differential control and so keep the cabin safe for rapid 
escape without risk of unconsciousness for a pilot from a 
too severe explosive decompression. This statement 
should apply for all flight conditions and not be limited 
to combat conditions alone. 


EMERGENCY PROCEDURES IN PRESSURE-CABIN AIR- 


CRAFT 


There are two critical points in the oxygen altitude 
scale that require warning for the air crewmen. The 
first is at 10,000 ft.—a warning to put on oxygen 
equipment. The second is at 40,000 ft., above which 
pressure breathing should be used. All tactical pres- 
surized-cabin aircraft use a cabin altitude warning 
switch to indicate these critical levels by means of a 
horn or pilot light. At the present time, all fighters and 
bombers of the Army Air Forces have the demand 
oxygen system. Those capable of flight above 
35,000 ft. are equipped with the pressure-demand sys- 
tem. 

A factor that must be given serious consideration 
while flying in pressure-cabin aircraft is the duration of 
useful consciousness for a pilot after decompression. 
Above 30,000 ft. the time of useful consciousness be- 
comes increasingly less. Skill and practice will there- 
fore be required by the pilot to connect his oxygen sys- 
tem. The allowable time for emergency action is illus- 
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trated in Fig. 3. On the basis of this chart it may be 
concluded that all crews of pressure-cabin tactical air- 
craft in combat areas should wear their oxygen masks 
for flights above 30,000 ft. 

Below 30,000 ft. the loss of flying efficiency as a result 
of a sudden decompression of the cabin is a matter of 
minutes rather than seconds. Curves indicating these 
tolerance limits are presented in Fig. 4. Below 25,000 
ft. the limits vary greatly with individual tolerance to 
anoxia and may be considered in the present chart as 
indicating a trend rather than a basis for design or 
operation. 
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based on data presented in reference 1. 


Use of cabin pressurization serves to reduce the re- 
quired oxygen capacity rather than to eliminate the 
need for the system. Cargo and passenger aircraft 
with pressure cabins, such as the Constellation (C-69) 
and Strato Cruiser (C-97), have continuous flow systems 
for both the passengers and the crews. 


It would be a fallacious policy for designers and en- 
gineers to plan, for postwar commercial use, cabin- 
pressure aircraft without oxygen systems, when the 
danger from anoxia to passengers and air crews at al- 
titudes above 10,000 ft. is realized. Although, for the 
crew, descent from 20,000 ft. to safe altitudes is some- 
times possible with safety, occasions may arise when 
weather conditions preclude sudden lowering of alti- 
tude. Under these conditions the danger to crew and 
passengers is serious. An important safety rule for all 
pressure-cabin aircraft is that supplementary oxygen 
must be available at all altitudes when flights above 
10,000 ft. are made. Present-type continuous flow 
oxygen systems, such as used in the Constellation (C- 
69), weigh only 225 Ibs. and cause a weight penalty of 
approximately 3/2 lbs. per passenger and crew. This 
weight penalty may be reduced by use of high-pressure 
supply oxygen or of liquid supply oxygen. A new type 
of emergency system that is showing considerable prom- 
ise experimentally is a chemical oxygen generator using 
superoxide or perchlorate. Walk-around oxygen bottles 
for stewardesses and crew members will also be neces- 
sary emergency items. Finally, the weight penalty 
of oxygen equipment for commercial aircraft as safety 
insurance is small compared with the pay load and cer- 
tainly small compared with the 1,200 Ibs. now required 
by the Air Transport Command on transoceanic flights 
for life raft and fire-protection equipment. 


Use of oxygen supplied by continuous flow systems 
below 18,000 ft. does not require any skill or indoctrina- 
tion of the passenger. Above 18,000 ft. it may be nec- 
essary to have an oxygen mask drill for passengers in 
pressure-cabin aircraft much the same as life-boat 
drills are held on shipboard. Emergencies, caused by 
loss of cabin pressure above 25,000 ft., will require 
quick donning of oxygen equipment by crew and pas- 
sengers, as the time reserve for useful consciousness is 
4 min. or less. An additional emergency consideration 
above 25,000 ft. is the possibility that aero-embolism 
may occur in some passenger or crew member. As a 
general practice, plane crews, including stewardesses, 
must be thoroughly indoctrinated in low-pressure cham- 
bers. 


At this point it should be noted that the current Civil 
Aeronautics Board Regulation (Amendment No. 41.23, 
dated 27 June, 1945, effective 1 August, 1945) requires 
an emergency oxygen supply for flight in pressure- 
cabin aircraft above 18,000 ft. only for the crew but not 
for the passengers. This regulation would, therefore, 
fall far short of the safety requirement for care of pas- 
sengers outlined in the previous sections. 
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MEDICAL ASPECTS OF CABIN 


COMFORT IN PRESSURE-CABIN AIRCRAFT 


At the present time there are few actual flight data 
on comfort requirements in pressure-cabin aircraft that 
give a real basis for a clear statement of design require- 
ments. In general, the factors that must be considered 
are: (a) temperature and humidity; (b) ventilation; 
(c) sound control; and (d) choice of isobaric. 

Fig. 5 illustrates the human comfort zone and the 
tolerance levels for hot and cold temperature for air 
crews and passengers dressed in ordinary everyday 
clothes. The comfort and tolerance levels are defined 
by the environmental temperature, which may be 
either the dry bulb temperature or the resultant tem- 
perature when radiation effects are considered, and the 
absolute vapor pressure of the environment. It will be 
noted that in the warm ranges humidity plays an im- 
portant role in comfort, while toward the cold range its 
effect is negligible. The air motion equivalent for this 
chart is 200 linear ft. per min. 

The ventilation requirement in a pressure-cabin air- 
craft must be sufficient, first, to maintain the necessary 
cabin temperature when either heating or cooling is 
required and, second, to remove odors, carbon dioxide, 
and noxious gases. There are three types of ventilation 
that may be used in pressure-cabin aircraft. The con- 
ventional type, a ram source directed into the cabin 
from the slipstream, is called the auxiliary ventilation 
system and would be used only in an emergency or 
when pressurization is not used. With pressurization, 
the air from the supercharger acts as the main ventila- 
tion source. In bombardment and fighter aircraft this 
source of ventilation may amount to 70 to 100 cu.ft. 
per person, which is more than adequate for the re- 
moval of obnoxious gases. For passenger aircraft the 
ventilated air from the superchargers can drop as low 
as 4 to 6 cu.ft. per person. Although the latter value is 
sufficient to remove carbon dioxide, it often is insuffi- 
cient to keep the odor and cigarette smoke level down. 
Under these circumstances, recirculation of air through 
air filters and deodorizers may be necessary for comfort. 
If the experiences of the railroads in conditioning their 
pullman and dining cars is of any value in aircraft, the 
total ventilation in aircraft should be at least 10 cu.ft. 
per person for nonsmoking sections and 20 cu.ft. per 
person for compartments where smoking is allowed. 
Deodorization is a specially important consideration if 
patients are carried for long periods in flight. How- 
ever, there are no quantitative data available to specify 
this requirement. For completeness, it should be 
noted that a part of the ventilated air from the super- 
chargers must be used for defrosting windows and pre- 
venting condensation on them. Since with large pas- 
senger loads dew points as high as 50°F. can be expected 
where ventilation is not used, proper insulating meth- 
ods must be carefully chosen. 

Sound control in pressurized-cabin aircraft is no more 
difficult a problem than for conventional aircraft and 
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A.C., Chief, Thermal Research Unit, Biophysics Branch, Aero 
Medical Laboratory, presents the relation between comfort and 
tolerance to heat and cold in terms of hours, environmental tem- 
peratures, and so forth. In this chart environmental tempera- 
ture is the resultant of air temperature, wall temperature, and 
effective solar radiant temperature. It is assumed (a) that per- 
sonnel are sitting and doing no more than light manual work, 
(b) that standard temperate climate clothing (one clo) is worn, 
and (c) that air motion is the equivalent of 200 ft. per min. linear 
velocity. The tolerance lines give the recommended maximum 
duration of exposure to the stated conditions. 


in some respects may be easier. Preliminary data to 
date indicate that pressurization tends to remove noise 
entirely from the high spectral region and confine the 
noise level to frequencies below 1,200 cycles. 

Aside from the effect of altitude on night and periph- 
eral vision noted previously, there are other consid- 
erations for maintaining isobaric cabins at the 5,000-ft. 
level or lower. Comfort may be defined as the state 
in which one is unconscious of adaptation to his envi- 
ronment. Discomfort, on the other hand, implies con- 
sciousness of the accommodative processes that are be- 
ing made by the body. Discomfort can extend also to 
a level of pain. From the temperature point of view, a 
person is comfortable when he feels neither hot nor cold. 
When he is hot or cold, he can accommodate physio- 
logically to his environment but still be uncomfortable. 
When acclimatization takes place, he is comfortable 
since he is no longer conscious of his accommodation. 
After he is exposed to the environment, permanent 
physiologic changes take place which are defined as 
adaptation. Human response to altitude is parallel to 
the response to heat and cold. Wher a human being, 
whose normal environment is sea level, is flown to 8,000 
ft., he must accommodate to this level by changes in 
respiration. He is conscious of an accommodative ef- 
fort. (This is a common experience of air passengers 
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from the east or west coasts when they get off for a 
stretch at Albuquerque.) Although he is not neces- 
sarily cyanotic, he has a sensation of air starvation when 
exercising. Passengers do not have a sense of comfort 
or complete freedom from this sensation above 5,000 
ft. until they are acclimatized. Normal healthy people 
can accommodate to altitudes as high as 12,000 ft. On 
the other hand, cyanosis has been observed at altitudes 
as low as 8,000 ft. in normal subjects acclimatized to sea 
level. 


During the present war there has been a great deal of 
transoceanic flying by the Air Transport Command, the 
Naval Air Transport Service, and tactical aircraft at 
altitudes of 8,000 to 9,000 ft. without regular use of 
oxygen. At first, the same air crews were used for the 
complete transoceanic run without change. Now, 
crews are rotated every 8 to 10 hours at intermediate 
stops and are given 24 hours of rest before completing 
another leg. This rotation practice has resulted in a 
great improvement in pilot efficiency. After 24 to 36 
hours of flying at 8,000 to 9,000 ft. with occasional stops, 
passengers often report extreme lassitude for several 
days on returning to their businesses. 


Although at 8,000 ft. a pilot can act without handicap 
in an emergency, long exposure of 8 to 12 hours leaves 
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him, as a rule, sleepy, unable to concentrate, and with 
a general lack of motivation. In reality, the 8,000-f 
isobaric for pressure-cabin aircraft was proposed back 
in 1934 when flights longer than 4 hours were not con. 
templated as a routine. The experience of this war at 
flight levels of 8,000 to 9,000 ft. for 8 to 12 hours has 
demonstrated that, for pressure-cabin aircraft, comfort 
of passengers and efficiency of crews demands the use, 
wherever possible, of isobarics of 5,000 ft. or lower. 

In summary, aviation is now long past the pioneering 
stage in which comfort of passengers and efficiency of 
crew were minor factors compared to performance and 
speed and altitude records. Aviation engineering has 
advanced so far that the price, in weight penalty and 
performance, for safety and comfort is small. This is 
especially true when it is realized that speeds of more 
than 300 m.p.h. and flights of 8 to 12 hours are common- 
place. Cabin pressurization is the first step taken by 
the aviation industry to include in the structural design 
of an aircraft a factor devoted entirely to the safety 
and efficiency of air crew and the comfort of the passen- 
gers. 
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The Stress-Area Method Applied to Frames 


GILBERT C. BEST* 
Consolidated Vultee Aircraft Corporation 


SUMMARY 


A procedure is worked out to facilitate the design of frames. 
The necessary alterations in frame section moments of inertia 
to obtain predetermined stresses at given points are evaluated. 
The principal saving of labor is due to the use of original frame 
sections—the added material being considered removed and 
replaced by loads causing identical stresses throughout the struc- 
ture. An example problem is given illustrating the method. 
An appendix gives additional -proofs to the relation between 
section changes and equivalent moments and an extension of the 
stress-area method applied to beams! to include the Jo’s of added 


areas. 


INTRODUCTION 


I THE DESIGN OF FRAMES many successive trials are 
usually required to obtain sections producing satis- 
factory stresses in extreme fibers. When any alteration 
is made, a redistribution of moments about the frame 
occurs which complicates the procedure. It is possible 
to simplify matters somewhat by considering—after 
an alteration has been made—that all added material 
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is removed, and, wherever such material was cut away, 
forces replace the stresses formerly existing across the 
cut. Stresses throughout the structure remain un- 
changed by this operation, and it is then possible to 
work with the original frame sections. The law of 
superposition of loads also conveniently applies with 
regard to the new loading. It is obvious that what- 
ever changes in stress occur are then due to such 
“stress-area” loads. It is also obvious that these 
loads will depend upon the added area and also upon 
the moment—and in some cases the axial load—at 
the section in question. 


Basic EQUATIONS 


The energy method will be used,’ and, as is customary 
in the analysis of frames, the deflections due to axial 
loads and shear loads will be neglected. For definite- 
ness assume a single cut. This gives the following 
equations: 


—" 


| (20) 
| 





in which X,,, Xp, X, are the moment, axial load, and 
shear at the cut, respectively. M,,, M,, M, are the 
moments caused by a unit moment, unit axial load, 
and unit shear load, respectively, applied at the cut. 
M, = static moment. 

It is evident that, if an alteration is made and equiva- 
lent moments replace added material, the only change 
will occur in the static moments M, in the right mem- 
bers of Eqs. (20). Also from consideration of the 
solution by determinants the above may be solved 
with only the changes in right members as constants 
and the resulting values—which will be termed AX,,, 
AX,, AX,—added to the original values for Xm, X,, 
X, to obtain the final solution. Or, if the moment 
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at some section before alterations are made equals M, 
M, = M, + aXm + OX, + cX, (21) 


where a, b, and c are constants depending upon the loca- 
tion of the section. 
After alterations are made the moment will be 
M = M,+ a(Xm t+ AXm) + 0(X, + AX5) + 
c(X, + AX,) 


or 


(22) 


M = M, + aAX,, + bDAX, + cAX, 


EQUIVALENT MOMENTS 


Reducing the frame problem to its essentials, it is 
sufficient to have given a line in a plane representing 
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the location of the neutral axis, a value of J at each 
point on this line, and a loading that is in equilibrium, 
to determine the distribution of moment about the 
frame. From the moment and the bending equation 
My/I the stress at a distance y from the neutral axis 
can then be determined. In this way a certain ‘‘stress 
field’’ exists about each point on the neutral axis. The 
problem in finding a moment equivalent in its effect to 
material which has been cut away is to find for each 
removal a moment to apply at the affected portion 
which will leave the stress field everywhere unchanged. 
In particular the bending stress at a distance y from the 
neutral axis at the point being altered must remain 
unchanged, and when this is the case stresses elsewhere 
are unaffected, since the equivalent moment replaces 
the effect of the removed material. 

Hence, letting R = Inew/Jorig. and m, = equivalent 
moment, 


if, = My RI orig. = (M + Ma) (y¥ Ty) (23) 
or 
M[(1/k) — 1] = m, 


or 


—[(k — 1)/k|M = m, (24) 


Two other proofs of this relation are given in the 
Appendix. 


APPROXIMATION OF k AND ™, 


Since it is necessary to know M before k and m, at 
a section can bé computed and m, in turn causes changes 
in the moment M, it is necessary either to use a process 
of successive approximations or to solve the problem 
algebraically. In general, an algebraic solution ap- 
pears unsatisfactory, since it becomes unduly com- 
plicated if more than one cut is involved. 

Given the moment and axial load at any section after 
alteration, the stress becomes 


Prot. = (My/kI) + [P/A(R)] (25) 


where P = axial load and A(k) = a function of k 
representing the area. Signs are to be determined in 
this formula by inspection. k is always positive by 
definition. 


SHIFT OF NEUTRAL AXIS 


Another effect of axial load occurs when a section is 
altered so as to cause its neutral axis to shift in the 
plane of the frame by an amount e. In this case an 
additional moment Pe exists at the altered section. 
This effect will also put moments upon adjacent sec- 
tions. All such moments should theoretically - be 
added to M in Eq. (24). Furthermore, to be precise, 
revision of the dimensions locating the neutral axis 
would be necessary. Usually it appears satisfactory, 
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however, to neglect such changes, since their effect is 
small. ; 

From Eq. (25) k can be determined for any selected 
stress ficot, once M and P are known. Then from 
Eq. (24) m, can be obtained. Substituting m, for 
M, in Eggs. (20) gives AX,,, AX,, AX,. Since the left 
members of the equations are unaltered, this does not 
involve extensive work. A closer approximation to 
M may then be obtained from Eq. (22) and the process 
repeated until convergence is obtained. Substitution 
of the kI’s so obtained for the original J’s at each 
altered section in the recomputation of frame moments 
produces the desired stresses in extreme fibers. An 
example problem follows. 


EXAMPLE PROBLEM 


The frame of Fig. 1 is loaded as shown. To compute 
static moments a cut is made at the bottom centerline. 
Positive moments will be defined as producing com- 
pression in the outside fibers. Positive axial loads 
produce tension. Since the loading is symmetrical, 
there will be no shear at the cut, and Eqs. (20) become 


Computation of Frame Moments 
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Computation of Frame Moments 
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TABLE 


15 


13 


11 
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M, = @) + 


Ff, 
— 62,480 
— 45,944 
— 17,568 
— 18,389 
—34,194 
—32,180 
— 8,840 
— 14,500 
— 6,649 
— 7,881 
—11,491 
— 14,630 
— 17,664 


@3) + Xn 


x10-% 


ds 


ds 


X10-3 yM, 
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— 15,095 
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— 85,000 
—60,000 


— 36,200 
— 17,500 


¥ 
80.0 


Sta. 





— 29,404 
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A; = 


A; => 
—2,018.6 


— 146,233 
X 108 


579,188 


10,467 


213.89 


X 108 











* Distance from neutral axis to extreme fiber. 


, [ds 5 ds _ ds 
xa f$ x, fy$ - fs 


+ (26) 


ds ds ds 
Eu ae iw pe a — | My 
J 7 x, > 7 fr 


since X, = 0, M, = 1, M, = 
and may therefore be omitted. 

All frame sections are considered to be symmetrical. 
Frame stations are measured in inches along the neutral 
axis from the top centerline. 

Computation of station moments and extreme fiber 
stresses is shown in Table 1. A word of explanation 
may be required regarding the integrating constants 
shown in Col. 2. These follow Simpson’s one-third 
rule from stations 0 through 100 and the trapezoidal 
rule thereafter. The constant at station 100 is the 
sum of both. This was necessary since the evenness 
of the spacing became interrupted at this station. 

It is desired to reduce the stresses at stations 0 and 
10 to 38,000 Ibs. per sq.in. by alteration of sections just 
inboard of station 20. To obtain this result sections 
will be both deepened and increased in moment of 
inertia. Letting the right members of Eqs. (26) be 
represented by A; and Ag, solutions for X,, and X, as 
linear functions of A; and A, are shown in Eqs. (27). 


—y, and E is constant 








Computation of Frame Movements 


213.89X, — 10,467X, = —A, 


ll 


10,467Xm — 579,188X > —A, 


10°X, = —40.4344, + 0.73071A2 l 
27) 
Xp = —708.5 Ibs. 
Xm = —25,234 in. Ibs. 


| 

| 

| 

(27 
| 10°X, = —0.73071A, + 0.014932A_ 7 

| 
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This is done to facilitate later work. If, as in the 
column analogy method, the origin of the x,y axes is 
chosen at the elastic center of the frame, the terms 
involving fy(ds/I) vanish and the solution is thereby 
simplified. This has not been done in this example, 
since it was desired to demonstrate a method having 
a more general application. 

The required section change constants k are com- 
puted in Table 2, it only being necessary to carry the 
approximation process through one extra cycle to 
obtain satisfactory convergence. Lines 7 and 8, etc., 
may be omitted in an actual computation. They were 
inserted here only to clarify the example. Table 3 
shows a recomputation of section moments and extreme 
fiber stresses, indicating that the desired results are 
obtained within the accuracy to which the computa- 
tions were carried out. 
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TABLE 2 
Computation of k’s 
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Item Ref. 
eo Col. 6, Table 1 1.5 1.12 
2 Znew 2.5 2.0 
ug Col. 4, Table 1 80.0 78.8 
4  ds/I Col. 7, Table 1 2.22 11.90 
5 . M, Col. 14, Table1 —53,554 —29,404 
| 6 k=@ xX @/@® xX 38,000 Eq. (25) 2.35 1.38 
| elt ihe DM —0.574 —0.275 
= 8 m=6)X@®@ Eq. (24) +30,740 +8,086 
&/ 9 — ma(ds/I) = @ x ® x @- 68.2 96.2 <A, = 164.4 AX, = +75 
2 |10  — may(ds/I) = (3) xX ® x 1073 5,456 7,581 Az = 13,037 AXm = +2,879 
= ju —y X AX, —6,000 —5,910 Ref. Eq. (27 
12 M=@+@+ dX, Eq. (22) —56,675  —32,435 
133 k=Q@ xX @/@ X 38,000 Eq. (25) 2.486 1.524 
xil4  —(k —1)/k —0.598 —0.344 
E15 ma = @ xX Eq. (24) +33,892 +11,158 
3/16 = ma(ds/I) = @ X @ X 10-8 +75.2 +132.8 A,=+2080 AX, = +944 
~ {117 may(ds/I) = (3) X G8) X 10-8 +6016  +10,465 A: = 16,481 AXn = +3,633 
&\18 ~-yxX AX, —7,552 —7,439 Ref. Eq. (27) 
19 M=G6+@®+ AX, Eq. (22 — 57,473 —33,210 
20 k=Q@ xX ®@/@® X 38,000 2.52 1.56 
21. Inew = @) X 3.78 1.75 
Appendix 1 or M the relation is general. Any frame can be cut into 
subframes that are capable of analysis by the column 
ADDITIONAL PROOFS OF EQUIVALENT MOMENT analogy method, loads and moments replacing the 
FORMULA stresses across the cuts. The above substitution of m, 
for material removed can then be made on each sub- 
Proof A frame without affecting moment distribution thereon, 
g 


Consider that a section is so altered as to cause its 
neutral axis to shift parallel to itself a distance h. 
Let Lo = Lorig. and I, = Iaiterea; Similarly for Ao. Then 
—dm, = fodA(y + h) = dA(y + h)(My/I,), where y 
refers to the neutral axis after alterations have been 
made and dA = an infinitesimal part of the added area. 

Then —m, = (M/I,)[/ydA + hf ydA ], where in- 
tegrations are to be performed over the added areas only. 

Then, since fy dA = hAy and I, = Ip + h?Ay + 
Sy? dA =I,+ Al = kIh, 

—m, = MAI/(Ih + Al) = M(k — 1)/k 


Proof B 


In the column analogy method any change in section 
corresponds to the addition or subtraction of equivalent 
area ords/EI, Any change in section will therefore re- 
sult in a change in static loading of M,(ds/EI) [(1/k) —1]. 
Removal of the added ‘‘area’’ working at ‘‘stress”’ 
M;, and replacing by an equivalent ‘‘force’’ as per the 
stress area method applied to beams! gives a change of 
M((ds/EI)((1/k) — 1]. Since M = M, +°M,, the 
change in load to be applied to the original section is 
M(ds/EI)[(1/k) — 1]. Proceeding by the method of 
this article—i.e., dealing with actual areas—the change 
in loading is given directly as m,(ds/EI). Comparing 
with the above shows m, = M[(1/k) — 1]. Since 
no restrictions were placed upon the section change 


and subframes can then be put back together into the 
original frame. 


Appendix 2 


EXTENSION OF THE STRESS-AREA METHOD APPLIED 
TO BEAMS TO INCLUDE J)’s OF ADDED AREAS! 


Since the slope of the stress distribution curve is 
constant because of the original assumption of simple 
bending, the slope will be the same for the small added 
area as for the whole section. Hence it is possible to 
write 

M,/I- a — Moz/Ivz (14) 
where Mo,, Jo, are the reacting moment and moment 
of inertia of the small added area about the x axis, and 
M,; = applied moment M, plus the moment due to the 
added area. The minus sign is due to the fact that 
Moz is a reaction and M, an applied moment. 

Therefore, when the loads due to the new area at 
point Py working at stress Fy are considered, the small 
moment Mo, must be added giving—at the c.g. of the 
original section— 


Meay = aoF ovo —_ Mor 
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Ao F oxo 


= —<doFo 
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TABLE 3. 
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TABLE 3. 
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Fy 
—38,142 


—38,174 
— 25,975 
— 9,708 
— 27,355 
— 27,456 
— 8,096 
— 11,264 
— 5,529 
— 6,635 
— 9,720 
—12,396 
— 14,974 


@) + Xn 


— 57,671 


—yX py 
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x<10-8 yM, = 


ue 
I 


M 


— 85,000 
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— 5,992 
— 36,027 
— 23,763 
— 35,250 


— 7,780 


— 6,732 
— 1,349 


— 74.9 
— 457.2 
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— 1,648 
x 108 
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X 10 
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— 21,391 in.lbs. 


—609 lbs. 


= +116,893,000. X, = 


10,022Xm — 543,993X, 


10,022X, = +1,648,000. 


208.27Xm — 


However Mo; = —M,Io,/Ip and Mz = Mz + Mo, + 
QoFoyo; therefore, letting Jo,/J, = wz gives: Mo, = 
[—tor/(1 + woz) ](Mz + aoFoyo) and the first of Eqs. 
(32) becomes 


loz lor 


Mx. = Foy — —_ WM, ee een Fivo = 
a0 Aol9Vo f+“ Yor AoloVo 
ao lor 
— Foyo - —— M, (16 
lt lt (16) 


Writing the «, in the numerator of the second term as 
lor = Toz/I- _ Boz(aohoz?/TI,) 


where /o, is the length of the added area in a direction 
perpendicular to the x axis and fo, is a dimensionless 
constant depending upon the shape of the added area— 


e.g., if rectangular So, = '/1.—gives 
ao M; | 
Mz, = ——| Fow — — ho,? 
oT -| oo — Boz 2 10 
ao M 17 
Mow ELF Oem] (OP 
Pap = —ayFo 


Substituting Eqs. (34) into the bending Eq. (1)! gives 
AFip = (—10'Fo + a10’)do (18) 


where 
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| k2BorMzhoz* ds Ri BoyM yhoy? | (19a) 
i+al, +a, l” 


For practical purposes the shape may be assumed 
rectangular and 1 + uw, and 1 + w, roughly approxi- 
mated or, if « is small, set equal to one. If set equal 
to one, Yio’ becomes yo. It should be noted that 
ayo’ does not necessarily equal ag’; nor does yi’ neces- 
sarily equal yo’, so symmetry is lost when the effect 
of the Jo’s is considered. 

In general, when the effect of the moments of inertia 
of added areas about their own centroids is considered, 
terms of the type (— 10’ Fo + a1w’)do replace terms of 
the type — y10/d in the general equations. 
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SIXTH INTERNATIONAL CONGRESS 
FOR APPLIED MECHANICS 


At a meeting initiated by Dr. J. C. Hunsaker and Dr. Th. von K4rman, acting as joint secretaries 
of the Fifth International Congress for Applied Mechanics, a committee was formed for the organiza- 
tion of the Sixth Congress. In line with the decision reached at Cambridge in 1938, it is proposed 
that the Sixth International Congress for Applied Mechanics be held in Paris, from September 22 to 
September 29, 1946. 

The invitations to the Congress are extended on behalf of 1’ Académie des Sciences de 1’Institut 
de France, la Direction des Relations culturelles, le Centre national de la Recherche scientifique, 
l'Institut de Mécanique de la Faculté des Sciences de Paris, la Société francaise des Mécaniciens, 
and |’Association technique Maritime et Aéronautique. 


The Congress will meet at the Sorbonne. 
The Congress will be divided into the following Sections: 


I. Structures. Elasticity. Plasticity. 
II. Hydro- and Aerodynamics. Hydraulics. 
III. Solid Dynamics. Vibration and Sound. Friction and Lubrication. 
IV. Thermodynamics. Heat Transfer. Combustion. Fundamentals of Nuclear Energy. 
Besides the papers presented in these Sections, a number of General Lectures will be given on 
subjects of current interest. The titles of these Lectures will be made known in a later notice. 


The French Organization Committee of the Congress is composed as follows: 


President: Henri VILLAT 
Secretary General: Maurice ROY 


Members: 
Emile BARRILLON Frédéric JOLIOT 
Louis de BROGLIE ; 
Albert CAQUOT Henri a 
Paul DUMANOIS Joseph PERES 
Jules DRACH Ernest VESSIOT 


Those who desire to become members of the Congress are requested to inform the Secretary 
General as soon-as possible of their intention to attend. They should also indicate at the same time 
whether they wish to present a paper, and in what Section. This is required in order to facilitate the 
preparation of the program and the issuing of further notices. 


Communications are to be addressed to the Secretary General of the Sixth International Con- 
gress for Applied Mechanics. 
Institut Henri-Poincaré, 11, rue Pierre-Curie, PARIS (V°) 
Checks or money orders can be made payable to the: 
Sixiéme Congrés International de Mécanique appliquée 
Institut Henri-Poincaré, 11, rue Pierre-Curie, PARIS (V°) 
It is recommended to make payments on the following account: Cheques postaux Paris N° 649 
(address given above). . 


The authors who wish to present papers are requested to submit the titles of their contributions 
as soon as possible. As at previous Congresses, contributions must be of scientific value and perti- 
nent to the program of the Congress. 


The authors are also asked to send to the Secretary General, before July 31, 1946, either the 
manuscript of their paper (maximum 5,000 words) or a summary of no more than 300 words. It is 
most desirable to have at least the summaries of all papers available to the members at the time of 
the Congress, in order to facilitate discussions. The Organization Committee hopes it will be able 
to have the summaries reproduced at the right time. However, it would be a great help if authors 
who have facilities to prepare a sufficient number of copies of their summaries, or manuscripts, would 
do so and send such copies to the Secretary or bring them to the meeting. 


The Proceedings of the Congress, containing the minutes of the sessions, the texts of the General 
Lectures, and the papers presented in the Sections, will be published with the least possible delay. 
The members of the Congress will have the right to buy these volumes by subscription. The surplus 
copies will be sold later on at a higher price. 


Adequate facilities will be provided for the members for their stay in Paris. Arrangements will 
be made known in due time. As at the previous Congresses, excursions and visits to various institu- 
tions, places of scientific, historic, and artistic interest will be arranged. 


The payment of an inscription fee of 300 francs is required from the members of the Congress; 
this fee can be paid immediately or not later than on the opening day of the Congress. Ladies ac- 
companying members will not be required to pay such fee, except those who wish to address the 
meetings or to submit papers. 
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eAn Appreciation 
To Corporate Members of the 


Institute of the Aeronautical Sczences 


Listed below are the leading aeronautical and industrial companies that make possible 
the publication of the Journal of the Aeronautical Sciences. 
hey enable the Institute to publish scientific and technical papers in the Journal and to 
devote its pages exclusively to editorial material. 
Dues om individual members which would normally be used to publish this scientific 
and engineering monthly are thus released for additional services to the profession and industry. 
The Aeronautical Engineering Review and the Aeronautical Engineering Catalog, together with 
the JourNAL, provide a comprehensive technical information service for the benefit of the air- 





craft and related industries. 





Academy of Aeronautics 
Aero Insurance Underwriters 
Acrojet Engineering Corporation 
Aeroquip Corporation 
Agawam Aircraft Products, Inc. 
— Company Division, Joyce Aviation, 
b. 
Aircraft Radio Corporation 
Aireon Manufacturing Corporation 
AiResearch Manufacturing Company 
AiResearch Manufacturing Company of 
Arizona, Inc. 
Airex Manufacturing Company, Inc. 
Allis-Chalmers Manufacturing Company 
Allison Division, General Motors Corporation 
Aluminum Company of America 
American Airlines System 
American Overseas Airlines, Inc. 
American Bosch Corporation 
American Phenolic Corporation 
Associated Aviation Underwriters 
—— pee & Manufacturers, Inc. 
Atlas Su Compan 
Baker Steel Ye Tit Company 
Beech Aircraft Corporation 
Bell Aircraft Corporation 
Bendix Aviation Corporation 
Bendix Products Division 
Bendix Radio Corporation 
Eclipse-Pioneer Division 
Export Division 
Friez Instrument Division 
Pacific Division 
Scintilla Magneto Company 
The BG Corporation 
Boeing Aircraft Company 
Boeing School of Aeronautics 
Breeze Corporations, Inc. 
Essex Tool & Die Company 
Federal Laboratories, Inc. 
Cal-Aero Technical Institute, A Division of Grand 
Central Airport Company 
Carrier Corporation 


Cessna Aircraft Company 
Chandler-Evans Corporation 
Chase National Bank of the City of New York 
The Cleveland Graphite Bronze Company 
The Cleveland Pneumatic Tool Company 
Automotive-Aircraft Division 
Clifford Manufacturing Company 
Consolidated Vultee Aircraft Corporation 
Nashville Division 
Stinson Aircraft Division 
Continental Motors Corporation 
Cox and Stevens Aircraft Corporation 
Culver Aircraft Corporation 
Curtiss-Wright Corporation 
Airplane Division 
Development Division 
Propeller Division 
Wright Aeronautical Corporation 
Doak Aircraft Company, Inc. 
Douglas Aircraft Company, Inc. 
E] Segundo Plant 
Long Beach Plant 
The Dow Chemical Company 
Dzus Fastener Company, Inc. 
Eastern Air Lines, Inc. 
Eaton Manufacturing Company 


Thomas A. Edison, Incorporated, Instrument ° 


Division 
Edo Aircraft Corporation 
The Electric Auto-Lite Company 
Electrol Incorporated 
Engineering and Research Corporation 
Ethyl Corporation 
The Fafnir Bearing Company 
Fairchild Camera & Instrument Corporation 
Fairchild Aerial Surveys, Inc. 
Fairchild Engine and Airplane Corporation 
Duramold Division 
Fairchild Aircraft Division 
Ranger Aircraft Engines Division 
Federal Telephone and Radio Corporation 
Firestone Aircraft Company 
Fleetwings Division, Kaiser Cargo, Inc. 
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Fletcher Aviation Corporation 
G & A Aircraft, Inc. 
General Aircraft Equipment, Inc. 
General Aviation Equipment Company, Inc. 
General Electric Company 
General Instrument Corporation 
General Motors Corporation 
AC Spark Plug Division 
Aeroproducts Division 
Buick Motor Division 
Cadillac Motor Car Division 
Chevrolet Motor Division 
Delco Products Division 
Delco-Remy Division 
Eastern Aircraft Division 
Fisher Body Division 
Frigidaire Division 
Harrison Radiator Division 
Research Laboratories Division 
Rochester Products Division 
The B. F. Goodrich Company 
The Goodyear Tire & Rubber Company 
The Gray Manufacturing Company 
Grumman Aircraft Engineering Corporation 
Guiberson Diesel Engine Company 
W. & L. E. Gurley 
Hanna Engineering Works 
Haskelite Manufacturing Corporation 
Hawaiian Airlines Limited 
The Hilliard Corporation 
The International Nickel Company 
Jack & Heintz, Inc. 
Jacobs Aircraft Engine Company 
Johns-Manville Sales Corporation 
Casey Jones School of Aeronautics, Inc. 
Kellett Aircraft Corporation 
Kenyon Instrument Company, Inc. 
Walter Kidde & Company, Inc. 
Kollsman Instrument Division, Square D Com- 
pany 
Langley Corporation 
Lavelle Aircraft Corporation 
Lear Incorporated 
Liberty Aircraft Products Corporation 
Link Aviation Devices, Inc. 
The Liquidometer Corporation 
Lockheed Aircraft Corporation 
Longines-Wittnauer Watch Company, Inc. 
Magnaflux Corporation 
The Marquette Metal Products Company 
The Glenn L. Martin Company 
The W. L. Maxson Corporation 
Warren McArthur Corporation 
McDonnell Aircraft Corporation 
Menasco Manufacturing Company 
Micromatic Hone Corporation 
Minneapolis-Honeywell Regulator Company 
Moore Drop Forging Company 
National City Bank of New York 
National Credit Office, Inc. 
The New York Air Brake Company 
Norma-Hoffmann Bearings Corporation 
North American Aviation, Inc. 
North American Aviation, Inc., of Texas 
Northrop Aircraft, Inc. 


Northwest Airlines, Inc. 
Otto Aviation Corporation 
Owens-Corning Fiberglas Corporation 
Pan American World Airways System 
The Parker Appliance Company 
Pennsylvania-Central Airlines Corporation 
Pesco Products Company Division, Borg-Warner 
Corporation 
Phillips Petroleum Company 
Pioneer Parachute Company, Inc. 
The Pure Oil Company 
Republic Aviation Corporation 
J. P. Riddle Company 
A. V. Roe Canada Limited 
John A. Roebling’s Sons Company 
Rohr Aircraft Corporation 
Roosevelt Field, Inc. 
The Ryan Aeronautical Company 
Sciaky Brothers 
Scott Aviation Corporation 
Shell Oil Company, Inc. 
Simmonds Aerocessories, Inc. 
Skydyne, Inc. 
Socony-Vacuum Oil Company 
Solar Aircraft Company 
Sperry Gyroscope Company, Inc. 
Square D are seme 
Standard Oil Company of California 
Standard Oil Company (Indiana) 
Standard Oil Company of New Jersey 
Swedlow Aecroplastics Corporation 
Teleflex Limited 
The Texas Company 
Thompson Products, Inc. 
Tinnerman Products, Inc. 
Titeflex, Inc. 
Transcontinental & Western Air, Inc. 
Triplett & Barton, Inc. 
Union Carbide and Carbon Corporation 
Bakelite Corporation 
Haynes Stellite Company 
Linde Air Products Company 
National Carbon Company 
United Aircraft Corporation 
Chance Vought Aircraft Division 
Hamilton Standard Propellers Division 
Pratt & Whitney Aircraft Division 
Pratt & Whitney Aircraft Corporation of 
Missouri 
Sikorsky Aircraft Division 
United Air Lines, Inc. 
United States Aviation Underwriters, Inc. 
United States Rubber Company 
Dominion Rubber Company, Ltd. 
The Variety Aircraft Corporation 
Vickers, Inc. 
Vidal Corporation 
The Waco Aircraft Company 
Warner Aircraft Corporation 
The Weatherhead Company 
Western Air Lines, Inc. ; 
Westinghouse Electric Corporation 
Weston Electrical Instrument Corporation 
Wyman-Gordon Company 
Young Radiator Company 























Services 
of the 


Aeronautical Archives 
of the 


Institute of the Aeronautical Sciences 





The services of the Aeronautical Archives are available 
to all members of the Institute, to Corporate Members, to 
advertisers in the AERONAUTICAL ENGINEERING RgeviEw 
and, under usual library limitations, to the public. Four 
specialized services are available. 


The Paul Kollsman Lending Library 


This lending library service makes available, without 
charge, the latest and more important aeronautical books. 
Members may request the loan of any aeronautical or 
technical book they wish to borrow. Through an ex- 
change agreement with the Engineering Societies Library, 
any book on general engineering may be borrowed from its 
" great collection of over 160,000 volumes. 
A photostating service is available at usual library rates. 
Applications for membership in the library and further 
information will be sent on request. 


The W. A. M. Burden Reference Library 


This reference library contains over 12,000 aeronautical 
books, magazines, pamphlets, and reports gathered from 
world-wide sources and is one of the most ccmplete 
acronautical libraries in the world. Material from this 
library is not available for loan but may be used for 
reference purposes at the Aeronautical Archives. 


The Pacific Aeronautical Library 


6715 Hollywood Boulevard 
Hollywood 28, California 


Established in cooperation with the aircraft companies 
the library serves. The leading aircraft companies in or 
near Los Angeles participate in its support and operation. 

This service library for aeronautical research is available 
to the public for reading privileges. Source material in- 
cludes aerodynamic and structural research reports, as well 


as books on drafting, production methods, history, and al- 
lied sciences. It furnishes books, periodicals, and pamphlet 
material to the participating aircraft companies to supple- 
ment their engineering libraries. 


Technical Information Service 


This service has experienced personnel under the super- 
vision of trained aeronautical engineers to compile any in- 
formation desired. The services range from listing special- 
ized reference books to the preparation of exhaustive 
bibliographies, digesting of reports, and general surveys of 
any acronautical subject. Some of the available services 
are: 


Bibliographies on any aeronautical subject. 

Reports on any aeronautical subject. 

Digests of aeronautical books, papers, periodicals, and refer- 
ences. 

Translations. 

Engineering investigations of special aeronautical subjects. 

Biographies of individuals engaged in aeronautics. 

Photostats of any aeronautical or general engineering mate- 
rial. 

Microfilms made on special order. 

Photographs made from the Institute’s photographic collec- 
tion. 

Drawings and tracings made. 


In addition to the services mentioned any commission 
which comes within the scope of the Service will be ac- 
cepted. Special arrangements may be made for work re- 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any material in the Aeronautical Ar- 
chives of the Institute may be ordered at standard photostat 
rates. 
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